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ABSTRACT 


^  Three  areas  of  explosive  testing  are  investigated 
in  this  thesis.  First^a  model  of  the  blast  wave  is 
described  and  a  minimum  range-bandwidth  product  is 
calculated  for  several  percentages  of  peak  over¬ 
pressure.  This  product  is  based  on  a  single-pole,  low- 
pass  filter  whose  cutoff  frequency  turns  out  to  be  a 
constant  divided  by  the  positive  impulse  duration  of  the 
ideal  blast  wave.  The  second  topic  of  the  thesis  covers 
three  ways  of  removing  the  noise  from  the  blast  waves. 

One  of  these  methods  is  to  use  analog  notch  filters  on 
the  data  to  remove  the  frequencies  that  produce  the 
noise.  These  notch  filters  take  a  considerable  amount  of 
time  to  set-up.  Therefore,  a  digital  filter  was  designed 
to  remove  any  frequency  component  whose  amplitude  is 
greater  than  the  corresponding  frequency  component  of  a 
comparable  ideal  blast  wave.  ^This  works  very  well,  but 
for  those  that  prefer  numerical^  methods  to  digital 


filtering,  a  least  squares  program  was  also  developed. 
This  "sliding  least  squares”  domputes  an  estimate  of  the 
data  based  on  the  fact  that  </ny  curve  can  be  represented 


by  a  series  of  straight  lines.  )  The  third  section 
outlines  specifications  for  a  digital  telemetry  system. 
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This  system  consists  of  a  Central  Control  Unit  (CCU)  and 
many  Field  Units  (FU) •  The  CCU  will  control  the  test 
from  detonation  to  the  final  collection  of  data. 
Basically  it  will  contain  a  mini  or  main^fcam^computer 
and  a  transceiver. ^The  FUs  (up  to  1000)  are  to  be  based 
on  a  microprocessor  and  each  FU  will  handle  up  to  eight 
transducers.  Each  FU  will  be  battery  powered  and  be 
capable  of  operating  entirely  on  its  own  during  the 
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ABSTRACT 


Three  areas  of  explosive  testing  are  investigated 
in  this  thesis.  First  a  model  of  the  blast  wave  is 
described  and  a  minimum  range-bandwidth  product  is 
calculated  for  several  percentages  of  peak  over¬ 
pressure.  This  product  is  based  on  a  single-pole,  low- 
pass  filter  whose  cutoff  frequency  turns  out  to  be  a 
constant  divided  by  the  positive  impulse  duration  of  the 
ideal  blast  wave.  The  second  topic  of  the  thesis  covers 
three  ways  of  removing  the  noise  from  the  blast  waves. 

One  of  these  methods  is  to  use  analog  notch  filters  on 
the  data  to  remove  the  frequencies  that  produce  the 
noise.  These  notch  filters  take  a  considerable  amount  of 
time  to  set-up.  Therefore,  a  digital  filter  was  designed 
to  remove  any  frequency  component  whose  amplitude  is 
greater  than  the  corresponding  frequency  component  of  a 
comparable  ideal  blast  wave.  This  works  very  well,  but 
for  those  that  prefer  numerical  methods  to  digital 
filtering,  a  least  squares  program  was  also  developed. 
This  "sliding  least  squares"  computes  an  estimate  of  the 
data  based  on  the  fact  that  any  curve  can  be  represented 
by  a  series  of  straight  lines.  The  third  section 
outlines  specifications  for  a  digital  telemetry  system. 


This  system  consists  of  a  Central  Control  Unit  (CCU)  and 
many  Field  Units  (FU) .  The  CCU  will  control  the  test 
from  detonation  to  the  final  collection  of  data. 
Basically  it  will  contain  a  mini  or  main  frame  computer 
and  a  transceiver.  The  FUs  (up  to  1000)  are  to  be  based 
on  a  microprocessor  and  each  FU  will  handle  up  to  eight 
transducers.  Each  FU  will  be  battery  powered  and  be 
capable  of  operating  entirely  on  its  own  during  the 
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PREFACE 


This  thesis  is  written  with  two  primary  purposes  in 
mind.  First  and  foremost,  it  is  to  fulfill  the  require** 
ments  for  the  Master  of  Science  degree.  Department  of 
Electrical  Engineering,  Dniversity  of  Utah. 

Second,  the  work  is  intended  as  a  guide  for  those 
who  must  conduct  the  explosive  tests  as  a  part  of  their 
regular  work.  The  use  of  the  material  in  this  thesis 
should  make  their  work  easier  and  more  routine. 

My  sincere  thanks  to  Keaton  Turner  for  his 
technical  advice  and  to  my  wife  Mary  for  her  patience 
and  understanding. 

May  1982  Theodore  J.  Moody 


INTRODUCTION 


Three  important  aspects  of  explosive  testing  are 
described  in  this  thesis.  First,  the  mathematical  model 
of  a  blast  wave  is  described  and  a  calculation  performed 
to  determine  the  minimum  frequency  response  for  the 
instrumentation  system  used  to  record  the  data.  Second, 
three  methods  are  discussed  for  filtering  the  data  to 
remove  the  unwanted  noise  added  to  the  signal  as  a 
result  of  the  blast.  One  method  is  strictly  analog, 
another  uses  digital  filtering  techniques,  and  the  third 
is  a  numerical  method  of  smoothing  the  data.  The  final 
section  of  the  thesis  describes  a  telemetry  system  for 
collecting  the  data. 


FREQUENCY  RESPONSE 


1.  Model  of  Blast  Have 

Two  primary  factors  determine  the  destructive  force 
of  an  explosion.  The  first  is  the  peak  magnitude  of  the 
blast  wave  and  the  second  is  the  duration  of  the 
positive  phase  of  the  blast  wave  (the  portion  of  the 
blast  wave  whose  magnitude  is  greater  than  the  ambient 
pressure).  A  blast  wave  can  be  defined  as  the  wave  of 
overpressure  which  radiates  from  an  explosion  and 
overpressure  is  the  guage  pressure,  above  ambient,  in 
the  blast  wave. 

The  primary  purpose  of  the  electronics  used  in  this 
type  testing  is  to  measure  the  peak  magnitude  and 
duration  of  the  blast  wave.  A  blast  wave  is  very 
similar  to  a  sound  wave,  except  it  travels  at  a 
supersonic  velocity,  while  the  sound  wave  travels  at  the 
local  speed  of  sound.  In  fact,  one  definition  of  a 
blast  wave  is  a  compression  of  air  traveling  at  a 
velocity  greater  than  the  local  speed  of  sound.  Once 
the  velocity  decreases  to  the  speed  of  sound,  it  is  no 
longer  a  blast  wave,  but  instead,  a  simple  sound  wave. 

Unlike  sound  waves,  the  leading  edge  of  a  blast 
wave  (shock  front)  is  discontinuous,  much  like  the  waves 
in  the  ocean.  Even  when  the  source  of  the  wave  is  not 
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discontinuous,  the  shock  front  will  be.  One  of  the 
properties  of  air  that  determines  the  speed  of  sound  is 
its  density;  the  more  dense  the  air  the  greater  the 
propagation  velocity  of  the  wave,  be  it  sound  or  blast 
wave.  Therefore,  the  blast  wave  travels  fastest  in  the 
most  dense  part  of  the  compressed  wave,  thus  forcing  the 
leading  edge  of  the  wave  to  be  discontinuous,  as  can  be 
seen  in  Figure  1.  As  the  blast  wave  passes  the 
transducer  it  appears  as  the  pressure/time  history  of 
Figure  2. 

It  appears  to  have  an  exponential  or  logarithmic 
decay  and  in  fact  can  be  modeled  by  the  following 
equation: 


P(t)  *  Pm  (1-t/to)  exp(-at/to) 

where  P(t)  is  the  instantaneous  overpressure  at  time  t. 
Pm  is  the  maximum  or  peak  overpressure  observed  when  t 
is  zero,  (to)  is  the  point  in  time  where  the 
overpressure  returns  to  ambient  pressure  and  (exp)  is 
the  base  of  the  natural  logarithms. 

The  decay  parameter  (a) ,  of  the  above  equation,  may 
be  regarded  as  an  adjustable  factor.  It  is  selected  to 
provide  the  proper  decay  rate  to  match  the  appropriate 
values  of  blast  impulse.  In  other  words,  it  adjusts  the 
slope  of  the  decay  between  t  »  0  and  t  ■  to.  For 


Direction  of  Travel 


figure  1.  Development  of  explosive  Blest  Wave.  The  explosion 
say  start  as  a  syasetrical  shock  front,  but  as  it  travels  the 
velocity  is  greater  in  the  region  of  higher  pressure.  There* 
fore,  the  wave  progresses  as  shown,  becoming  discontinuous. 
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Figure  2.  Ideal  Blast  Wave.  The  curves  A  and  B  show  the 
effect  of  changing  the  decay  parameter  (a)  In  the  Frielander 
Approximation. 


composition  B  explosives  (the  type  in  the  AGM-65 
missiles),  using  a  ■  1.0  gives  a  very  close  fit  to  the 
experimental  data.  Therefore,  for  the  remainder  of 
this  paper  a  •  1.0  will  be  used.  This  equation1  is 
commonly  called  the  Prielander  Approximation. 

To  use  this  model,  one  must  first  know  the  peak 
overpressure  of  the  signal  and  the  time  for  return  to 
ambient  pressure.  To  determine  these  parameters,  one  can 
use  the  equations  developed  by  the  Department  of  Defense 
Explosive  Safety  Board  (DODESB)  which  are  based  on  an 
extensive  amount  of  experimential  data  collected  over 
years  of  testing.  These  equations*  are: 

for  0.5  <  y  <  440: 

In  (Pm)  «  7.0452041  -  1.6277561(ln(y) )  -  0.27399088* 
(ln(y))2  -  0 ,065773136* (In (y) ) 3  + 

0 .0065412563* (In (y) ) *  +  0.048236359* 

(ln(y) ) 5  -  0 .020072553* (In (y) ) ®  + 

0 .0030190449* (In (y) ) 7  -  0.00015984026* 
(ln(y) ) 8 

for  40  <  y  <  1000: 

Pm  -  (226.61762)/y1*4065913 
y  ■  range  /  (explosive  weight) 1//3 

1G.F.  Kinney,  EXBlSalVft.  fihafiJLB-JiLAi£r  (The 
MacMillan  Company,  New  York,  1962),  pp.  (3-87 . 

2 

L.  Giglio-Tos  and  T.E.  Linnenbrink,  "Air  Blast 
Pressure  Measurement  Systems  and  Techniques,”  Minutes. 
15th  Department  of  Defense  Explosive  Safety  Board 
Seminar,  September  1973,  p.  1400. 
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if  y  <  2.048: 

Impulse  ■  range  *  antilog (g) 

where: 

g  -  2.4177* (log (y) ) 2  -  1.5278*(log(y) )  +  1.30103 
if  y  >  2.048: 

Impulse  ■  range  *  42.5733  /  (y1*®29) . 


Since  the  impulse  is  actually  the  integral  of  PCt) 
from  t  ■  0  to  t  ■  to,  one  can  analytically  solve  the 
integral  for  (to)  with  the  following  result: 

to  •  0.001  *  impulse  *  2.718281828  /  Pm. 

The  result  of  these  equations  can  be  seen  in  the 
theoretical  plot  of  P(t)  for  a  TNT  explosive  weight  of 
100  pounds  at  a  range  of  25  feet.  Figure  3.  As  the  plot 
shows,  p(t)  actually  has  a  negative  component,  which 
should  not  be  surprising  due  to  conservation  of  matter 
(A  compression  of  air  followed  by  a  rarefaction  of  the 
air) . 

Since  the  blast  wave  is  discontinuous  at  the 
leading  edge,  one  would  expect  the  Fourier  transform  to 
have  an  infinite  spectrum,  which  is  true.  However,  in 
the  next  section  it  will  be  shown  that  with  the  proper 
use  of  filters,  one  can  predict  the  actual  peak 


amplitude  of  the  blast  wave  given  the  filter  response 
and  the  frequency  cutoff  of  the  low-pass  filter. 

2.  Frequency  Requirements 

3 

When  the  DODESB  wrote  their  directive  for  explosive 
testing,  they  relied  on  the  work4  of  L.  Giglio-Tos  and 
T.  E.  Linnenbrink  for  frequency  requirements  of  the 
data  acquisition  system.  In  developing  the  frequency 
requirements.  Hr.  Giglio-Tos  and  Hr.  Linnenbrink  used  a 
very  simple  model  for  the  blast  wave: 

P(t)  ■  Pm  *  (exp(-ct)) 

which  explains  why  their  curves.  Figure  4,  do  not  have  a 
direct  dependency  on  range.  This  dependency  comes  only 
when  one  takes  into  account  the  duration  of  the  positive 
portion  of  the  blast  wave,  (to),  as  described  in  the 
previous  section.  They  explain  that  their  curves  are 
based  on  an  explosive  weight  of  one  pound  and  to 
extrapolate  the  curves  to  other  explosive  weights  one 
should  divide  the  bandwidth  required  for  one  pound  by 
the  cube  root  of  the  actual  charge  weight.  This  is  an 

^Department  of  the  Air  Force  Technical  Order  TO  11A- 
1*47,  Department  of  Defense  Explosive  Hazard  Classifi¬ 
cation  Procedures.  Harch  1981,  para.  6-2. C. 

4L.  Giglio-Tos  and  T.E.  Linnenbrink,  "Air  Blast 
Pressure  Heasurement  Systems  and  Techniques,"  Minutes. 
15th  Department  of  Defense  Explosive  Safety  Board 
Seminar,  September  1973,  pp.  1359-1402. 
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R 

unproven  extension  of  the  scaling  law  which  has  been 
proven  by  others.  However,  OODESB  recommended  a  mini¬ 
mum  system  frequency  response^  of  20  kHz  for  all  tests. 

A  much  more  accurate  bandwidth  requirement  can  be 
obtained  by  starting  with  the  Frielander  approximation 
instead  of  the  simple  exponential. 

The  Laplace  transform  of  this  approximation  is: 
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Using  this  procedure  gives: 

Y(8)  -  P(s)H(s) 

Substituting  equations  (1)  and  (2)  for  PCs)  and  H(s): 

Xial  „  _ a _  _  _ a/ to _ , 

Pm  (s  +  a)  (s  +  1/to)  (s  +  a)  (s  +  1/to)  . 

Expanding  these  two  terms  with  Partial  Fractions  yields: 

xifii  .  ..ato  j _ i _ i_\ 

Pm  ato-1  Is  +  1/to  s  +  a) 


and  transforming  back  to  the  time  domain: 


y (t)  »  ato  lexpC-t/to)  -  exp(-at)! 

Pm  ato-1 


ato  -  [exp (-at)  -  exp(-t/to)l 
<ato-l) 1 


ato  It/tol [exp(-t/to) J  • 
ato-1 

Rearranging  and  combining  terms  yields  the  data  as  would 
be  recorded  by  the  instrumentation  system: 
y(t)  ■  Pm  Uato/ (ato-1)  1 1  (1-t/to)  (exp(-t/to) )  -exp(-at)  1 

+  tato/(ato-l)2] [exp(-t/to)-  exp(-at)] ).  (3) 

As  a  simple  check  of  this  result  let  (a)  approach 
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infinity  in  equation  (3).  This  yields  the  original 
Frielander  approximation,  which  it  should,  since  letting 
the  low-pass  filter  have  an  infinite  cutoff  frequency  is 
the  same  as  B(s)  ■  1,  or  not  having  a  filter  at  all. 

The  objective  of  this  mathematical  exercise  is  to 
determine  the  cutoff  frequency  of  the  filter  which  will 
produce  a  specific  percentage  of  y(t).  Therefore,  the 
next  logical  step  should  be  to  take  the  derivative  of 
y(t)  with  respect  to  (t).  Then  set  the  derivative  equal 
to  zero  and  solve  for  (t) .  This  would  determine  the 
time  (t)  that  the  maximum  occures.  The  next  step  would 
be  to  substitute  this  value  of  (t)  into  equation  (3)  and 
solve  for  the  value  of  (a)  that  produces  the  desired 
percentage.  However,  equation  (3)  is  transidential  and 
as  such  cannot  be  solved  by  normal  means. 

ne  method  that  works  fairly  well  is  to  guess  a 
value  for  (a)  and  increment  (t)  until  a  maximum  is 
reached.  If  it  is  not  the  desired  maximum,  then  a  new 
value  of  (a)  is  tried  and  the  process  is  repeated  until 
the  desired  percentage  is  achieved. 

Through  trial  and  error  calculations,  I  found  the 
cutoff  frequency  in  hertz,  (a),  to  have  the  relation: 

a  ■  k/to 


where  (k)  varies  according  to  the  desired  percentage  of 
Pm.  This  variation  is  shown  in  the  following  table. 
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Percentages  and  Corresponding  Constants 

Percentage 

k 

99 

205.435 

96 

89.290 

97 

54.1365 

96 

37.6695 

95 

26.2816 

90 

11.1493 

The  program  Response  (Appendix  A)  was  an  aid  in 
performing  this  iteration.  It  asks  the  operator  for  a 
pole  multiplier  (pmult) ,  a  constant  equal  to  twice  the 
value  of  (k) .  Then  the  program  selects  a  value  for 
lambda  (range  divided  by  the  cube  root  of  the  explosive 
weight)  and  computes  the  necessary  values  of  impulse. 

Pm,  and  (to) . 

With  (to)  computed,  the  radian  value  of  the  cutoff 
frequency  is  computed  by: 

a  ■  (pmult)* (pi) /to. 

Then  the  time  interval  between  t  •  0  and  t  »  to  is  divid¬ 
ed  into  10  000  increments  and  equation  (3)  is  tested 
with  each  increment  starting  with  t  ■  0  and  stopping 
when  a  maximum  is  found. 
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In  the  calculation  of  y(t) ,  Pm  ■  1.0,  therefore, 
the  maximum  value  of  y(t)  is  actually  the  percentage  of 
maximum  peak  overpressure  of  the  ideal  blast  wave.  This 
value  of  y(t)  along  with  the  corresponding  values  of 
bandwidth  (in  hertz)  and  lambda  are  output  to  files  for 
future  reference. 

The  program  then  returns  to  select  a  new  value  for 
lambda  and  the  process  starts  over  again.  This  loop 
ends  up  computing  200  increments  of  lambda  and  the 
corresponding  values  of  y(t)  and  bandwidth. 

The  program  assumes  a  value  of  range  equal  to  one 
foot.  Therefore,  the  bandwidth  that  is  computed  by  the 
program  is  actually  a  range-bandwidth  product.  The 
range-bandwidth  product  is  plotted  in  Figure  5  against 
peak  overpressure  for  90  and  95  to  99  percent  of  peak 
overpressure.  Comparing  Figures  4  and  5,  over  the  range 
of  1-1000  psi,  a  large  discrepancy  can  be  seen.  This 
difference  is  due  primarily  to  the  simplistic  model  that 
was  used  to  develope  the  curves  of  Figure  4. 

There  appears  to  be  a  good  correlation  between  the 
curve  published  by  DODESB  for  peak  overpressure  vs. 
scaled  distance  (lambda)  and  the  equations  that  were 
derived  from  this  curve.  The  DODESB  curve  is  shown  in 
Figure  6  (P  curve)  and  the  curve  plotted  from  the 
equations  is  shown  in  Figure  7  (P  curve) • 

The  equations  I  developed  for  the  impulse,  appear 
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DODESB  Overpressure  and  Impulse  vs.  Distance  From 
TNT  Hemishpere 


Figure  7.  Calculated  Overpressure  and  Impulse  vs.  Distance 
From  TNT  Hemisphere 


to  match  the  curve  of  Figure  6  quite  well  over  the 
entire  interval.  These  equations  were  required  for  use 
in  the  computer  programs. 

The  range-bandwidth  product  is  plotted  vs.  lambda 
for  90  and  95  to  99  percent  of  peak  overpressure  in 
Figure  8.  This  plot  is  possibly  more  useful  to  the 
engineer  than  Figure  5  because  other  parameters  he  has 
to  calculate  are  also  in  terms  of  lambda. 

In  determining  the  required  bandwidth  or  frequency 
response  for  a  particular  transducer ,  one  must  first 
estimate  lambda,  read  the  range-bandwidth  product  from 
Appendix  F  ,  and  then  divide  this  product  by  the  distance 
the  transducer  will  be  located  from  the  center  of  the 
explosion,  ground  zero.  If  a  single-pole  filter  is  used 
for  each  transducer,  the  actual  peak  overpressure  could 
be  calculated. 

As  an  example,  assume  that  a  transducer  is  to  be 
placed  25  feet  from  a  100  pound  (TNT  equivalent) 
explosion.  Lambda  would  equal  25  divided  by  the  cube 
root  of  100  or: 

lambda  -  25/4.6416  -  5.3861. 

Then  from  Appendix  F,  for  95  percent,  the  range- 
bandwidth  product  would  be  188  440.  Dividing  this  by  25 
feet  produces  a  required  bandwidth  of  7538  Bz. 


Range-Bandwidth  Product  vs.  Lambda 


The  filter  for  this  transducer  would  be  set  for  a 
cutoff  frequency  of  7540  Bz. 

If  after  the  test,  it  was  found  that  this 
transducer  measured  a  peak  overpressure  of  34  psi,  then 
the  actual  peak  overpressure  could  be  calculated  by 
dividing  the  34  psi  by  0.95,  ort 

Pm  ■  Ymax  /  % 

-  34  /  0.95 
*  35.8  psi. 

Note  also,  that  for  a  one  pound  charge,  the  instru¬ 
mentation  system  with  a  bandwidth  of  20  kHz  will  record 
at  least  95%  of  the  peak  overpressure  if  the  transducer 
is  located  eight  feet  from  ground  zero  and  99%  only  if 
located  over  200  feet  away.  Likewise,  if  the  charge  is 
1000  pounds  then  the  minimum  range  for  the  transducer 
would  be  18  feet  to  assure  recording  at  least  95%  of 
peak  overpressure. 

The  above  makes  two  obvious  assumptions,  first  that 
the  TNT  equivalent  explosive  weight  of  the  charge  being 
tested  can  be  accurately  estimated.  If  this  were  the 
case,  there  would  be  little  need  for  the  test  in  the 
first  place.  However,  for  lambda  in  the  range  of  2-60, 
the  curves  of  Figure  8  are  relatively  constant  and  a 
poor  guess  of  explosive  weight  would  make  less  than  1% 
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change  in  the  amplitude  of  the  recorded  peak.  For  this 
reason,  the  transducers  should  be  located  in  this 
region,  when  possible. 

The  second  assumption  is  that  the  received  signal, 
y(t),  is  free  of  noise  and  has  an  ideal  shape.  This  is 
very  unrealistic  assumption  and  dealing  with  the  noise 
is  a  major  part  of  reducing  the  data  to  useable 
information.  It  is  the  subject  of  the  following 
section. 


FILTERING  METHODS 


1.  Sources  of  Noise 

There  are  many  types  of  noise  (defined  here  as 
unwanted  signal)  added  to  the  blast  wave  by  the  time  it 
is  recorded.  Among  them  are  the  noises  commonly 
encountered  by  an  instrumentation  system  such  as  white 
noise  and  1/f  noise.  These  can  be  dealt  with  by 
standard  means. 

However,  there  are  noises  that  are  peculiar  to 
explosive  testing.  These  noises  create  a  real  problem 
for  the  engineer  trying  to  produce  reliable  data  for  the 
agency  requesting  the  test. 

One  of  these  noises  contaminating  the  blast  waves 
is  called  Electro-Magnetic  Pulses  (EMP) ,  caused  by  the 
detonation  of  the  explosive.  Any  time  there  is  an 
explosion,  energy  is  released  at  all  frequencies  from  dc 
to  those  in  the  spectrum  of  light.  A  certain  percentage 
are  in  the  radio  frequency  range  and  are  coupled  into 
the  instrumentation  system  through  the  wiring. 

Preventing  the  EMP  from  being  recorded  is  extremely 
difficult  and  expensive.  However,  the  EMP  is  relatively 
easy  to  seperate  from  the  data  since  it  only  occurs  at 
the  instant  of  the  explosion  and  the  blast  wave  takes  a 


23 


finite  amount  of  time,  usually  several  milliseconds,  to 
reach  the  transducers.  Therefore,  the  EMP  noise  can  be 
literally  erased  or  ignored.  This  EMP  noise  pulse  can 
even  be  helpful  when  used  as  a  start  time  to  measure  the 
time  of  travel  of  the  blast  wave  between  ground  zero  and 
the  transducer. 

Another  form  of  noise  is  very  large  spikes  with  a 
very  short  duration.  Since  the  signal  from  the 
transducer  is  zero  99%  of  the  time  and  only  takes  a  non¬ 
zero  value  during  the  passage  of  the  blast  wave,  the 
average  output  signal  power  of  the  transducer  is  very 
small.  At  the  same  time  noise  (white)  has  an  average 
amplitude  that  may  be  as  high  as  150  millivolts.  Mr. 
Schwartz  has  shown1  that  the  output  noise  power  is  not 
only  proportional  to  the  average  input  amplitude,  but 
also  the  power  is  directly  proportional  to  the  bandwidth 
of  the  system.  Therefore,  the  signal  to  noise  ratio 
(SNR)  for  this  type  of  system  is  going  to  be  inversely 
proportional  to  the  bandwidth  and  quite  small  to  start 
with. 

One  of  the  advantages  of  using  an  FM  recording 
system  is  the  ability  to  record  very  low  frequencies 
accurately.  However,  below  an  SNR  of  about  13  dB  the 
noise  "quieting*  properties  start  to  deteriorate  rapidly 

1M.  Schwartz,  Information  Transmission.  Modulation 
and  Noise.  3rd  ed.,  (McGraw-Hill  Book  Company,  New  York, 
1980) ,  p.  351. 
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and  below  8  dB  an  AM  system  actually  has  better  noise 
characteristics.  At  even  lower  values  of  SNR,  the 
analysis  falls  appart  completely  because  it  is  based  on 
a  large  SNR. 

Since  the  signal  is  small  to  begin  with,  increasing 
the  bandwidth  decreases  the  SNR  to  the  point  that  the 
system  will  see  noise  that  would  not  be  there 
otherwise.  Theoretically  the  noise  could  dominate  the 
system. 

I  recommend  the  use  of  filters  on  each  of  the 
transducers  with  the  lowest  possible  cutoff  frequency  to 
reduce  this  type  of  noise.  Another  way  to  reduce  the 
noise  would  be  to  not  use  the  FN  system  at  all,  instead 
one  could  use  an  AM  system.  Even  better  would  be  a 
digital  system  like  the  one  described  in  the  last 
section  of  this  thesis. 

The  most  significant  problems  of  noise  occur  when  a 
blast  fragment  passes  near  a  transducer.  These  fragments 
cause  a  sonic  boom  type  shock  wave  to  be  superimposed  on 
the  blast  wave.  The  result  is  a  signal  in  which  one  can¬ 
not  determine  the  peak  of  the  blast  wave  nor  the  positive 
phase  impulse.  See  Figure  9  for  an  example. 

Analyzing  and  removing  this  type  of  noise  is  the 
subject  of  the  following  sections. 


2.  Spectrum  Analysis 

Before  attempting  to  filter  the  data,  I  wrote  a 
Fortran  program  to  determine  the  power  spectrum  of  a 
signal.  This  program  (Appendix  B) ,  Spectrum,  reads  the 
signal  from  a  file  and  computes  the  magnitude  and  power 
spectrum  of  the  Fourier  transform.  Spectrum  then  files 
each  in  a  seperate  file  for  later  plotting. 

Figure  10  shows  the  power  spectrum  of  the  actual 
test  data  in  Figure  9,  and  Figure  11  shows  the  power 
spectrum  of  a  comparable  ideal  blast  wave.  Comparing 
the  two,  we  can  see  roughly  what  would  have  to  be  done 
to  make  the  actual  data  have  the  same  shape  aB  the 
ideal. 

Many  of  the  frequencies  would  have  to  be  reduced 
and  many  others  should  be  increased.  One  cannot 
totally  remove  the  frequencies  that  exceed  the  ideal 
power  spectrum  (comb  filter) ,  because  a  certain 
amplitude  of  those  frequencies  is  required  to  make  up 
the  ideal  wave.  Likewise,  how  much  should  the 
frequencies  below  the  ideal  be  increased?  if  one  makes 
the  power  spectrum  of  the  actual  data  match  the  ideal 
exactly,  then  nothing  has  been  gained.  However,  if  the 
actual  data  spectrum  could  be  made  to  have  the  same 
general  shape  as  the  ideal  with  the  amplitude 
proportional  to  the  actual,  then  the  data  would  appear 
identical  in  shape  to  an  ideal  blast  wave. 
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Each  transducer  and  each  test  produces  a  different 
set  of  noise  frequencies  to  be  removed  from  the  data. 
Therefore,  designing  a  filter  to  remove  the  noise  before 
it  is  recorded  would  be  very  difficult,  perhaps  an 
adaptive  filter  could  be  used.  A  simpler  method, 
though,  might  be  to  record  the  noisy  data  and  then  try 
to  remove  the  noise, 

3.  Analog  Filters 

The  data  I  had  to  work  with  was  stored  on  analog 
tape,  recorded  at  120  inches  per  second  (ips).  This 
tape  recorder  will  play  back  data  at  any  of  several 
slower  speeds.  A  playback  speed  of  15  ips  was  chosen, 
this  speed  reduces  all  frequencies  by  a  factor  of  eight 
and  increases  all  times  by  the  same  factor  of  eight. 

Figure  12  is  the  filter  circuit  I  designed  to  help 
remove  the  sonic  boom  noise  from  the  blast  wave.  It 
incorporates  a  single-pole  filter  because  the  data  I  was 
working  with  had  not  been  filtered  previously. 

The  section  of  the  circuit  consisting  of  Rl,  R2, 
and  Cl  through  C4  is  the  low-pass  filter  and  has  a 
variable  cutoff  frequency  of  125  to  10  000  Hz  as 
indicated  in  the  figure.  IC1  is  a  voltage  follower  to 
isolate  the  low-  pass  filter  from  the  notch  filter, 
which  is  the  remainder  of  the  circuit. 

A  standard  active,  band-pass  filter  has  a  transfer 
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SWITCH  1,0.  Fc  (Hi)  R1*R2  (k) 

1  125-500  127-32 

2  625-1000  51-32  (Center  Frequency  of  notch  1g 

3  1250-5000  127-3 2  variable  over  tbe  Range  of 

4  6250-10,000  51-32  25-250  Hz) 

Adjustaente: 

1.  Connect  sine  wave  generator  to  Via  and  scope  to  Vout. 

2.  Set  R4  to  aaxlaua  real* tenet. 

3.  Set  elne  wave  generator  to  deaired  low-paee  cutoff  frequency. 

4.  Select  proper  capacitor  with  awitch. 

5.  Adjust  Rl  until  output  is  0.707  tiaes  tbe  input. 

6.  Set  cHne  wave  generator  to  desired  notch  frequency. 

7.  Adjust  R4  for  ainiaua  output. 

8.  Adjust  RJ  for  dsslred  attenuation. 


Figure  12.  Adjustable  Filter  Schematic 


function: 


H (s)  -  _ n 

8  +  a  (Wo)  8  +  (wo)  . 

Mote  the  negative  sign.  Adding  the  orignal  signal  to  the 
output  of  the  band-pass  filter  makes  the  circuit  a  notch 
filter  with  the  transfer  function: 

h(s)  «  l  - _ 5 _ .LflalalttglB--, 

s*  +  a  (Wo)  s  +  (wo)  . 

The  reason  for  using  this  arrangement  is  that  the  midband 
voltage  gain  and  the  center  frequency  of  the  notch  are 
easily  adjustable.  R3  adjusts  the  midband  gain  and  R4 
adjusts  the  center  frequency.  The  center  frequency  can 
be  adjusted  over  a  full  decade,  25  to  250  Hzr  while 
maintaining  a  constant  notch  width  of  10  Hz.  Decreasing 
C5  and  C6  by  an  order  of  magnitude  will  increase  the 
center  frequency  and  the  notch  width  by  an  order  of 
magnitude.  The  bode  plots  of  Figures  13  through  16  show 
the  versitility  of  this  circuit. 

Figure  17  shows  the  results  of  filtering  the  actual 
test  data  with  the  notch  set  at  38.75  Hz  and  the  low- 
pass  cutoff  set  at  149.7  Hz.  Recall  the  data  is  being 
played  back  at  a  reduced  speed  so  this  would  correspond 
to  having  a  notch  at  310  Hz  and  a  low-pass  cutoff  of 
just  under  1200  Hz  in  real  time. 

Even  though  this  notch  removed  the  major  noise 
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Figure  17.  Data  After  Analog  Filter  (Notch  ■  38.75  Hzf 
Low-Pass  ■  149.7  Hz).  Top  trace  is  Raw  Data  and  bottom  trace 
is  filtered  data.  The  horizontal  direction  is  timer  approxi¬ 
mately  1.88  msec,  per  small  division.  Vertical  Scale  is  volt¬ 
age  each  small  division  is  approximately  80.3  mV.. 
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signal  from  the  data,  additional  notch  filters  would  be 
required  to  smooth  the  data  properly.  These  would  be 
set  for  the  second  and  third  harmonics  of  the  original 
notch.  Also  a  60  Hz  (480  Hz  real  time)  notch  seemed  to 
remove  some  of  the  noise.  In  all  then  one  would  need  a 
minimum  of  four  notch  filters  in  addition  to  the  low- 
pass  filter,  to  properly  condition  this  signal. 

The  set-up  and  adjustment  of  this  many  filters  would 
be  very  time  consuming.  I  estimated  that  with  practice 
and  the  aid  of  a  spectrum  analyzer,  a  technician  could 
process  one  signal  in  just  under  under  four  hours. 
Realizing  that  a  given  test  may  produce  20  to  50  such 
signals  shows  that  the  man-hour  cost  would  probably  be 
prohibitive.  Attempting  to  do  this  kind  of  filtering 
without  the  aid  of  a  spectrum  analyzer  would  be 
unthinkable. 

A  much  more  accurate  and  faster  method  of  filtering 
the  data  would  be  to  digitize  the  data  and  then  use  the 
computer  to  remove  the  noise. 

Although  numerical  smoothing  is  a  form  of  digital 
filtering,  I  have  separated  them  in  the  following 
sections.  In  this  thesis,  digital  filtering  indicates  a 
frequency  domain  manipulation,  while  numerical 
smoothing  operates  only  on  the  time  domain  Bignal. 
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4.  Digital  Filter 

Fourier  analysis  has  been  around  £or  a  long  tine, 
but  it  has  only  been  since  the  advent  of  the  large 
memory,  high  speed,  digital  computers  that  the  use  of 
this  method  for  filtering  has  been  feasible. 

The  Fast  Fourier  Transform  (FFT)  is  an  efficient 
algorithm  for  computing  the  discrete  Fourier  transform. 
Using  the  FFT,  one  can  compute  the  Fourier  transform  of 
a  signal,  decide  what  frequency  components  in  the  signal 
are  not  wanted,  remove  all  or  part  of  them,  and  then 
return  the  signal  to  the  time  domain  as  a  filtered 
version  of  the  original  signal. 

Although  a  simplified  description  of  digital 
filtering,  this  is  in  essence  what  a  digital  filter 
accomplishes.  The  big  advantage  of  a  digital  filter  over 
its  analog  counterpart  is  the  ability  to  rapidly  and 
precisely  change  the  characteristics  of  the  filter  to 
match  the  situation.  An  analog  filter  has  to  make  a 
change  in  hardware  while  the  digital  filter  only  has  to 
change  some  bits  in  memory. 

For  the  particular  problem  of  filtering  the  blast 
wave  data,  the  first  task  was  to  convert  the  analog 
signal  to  a  digital  signal  (Figure  18) .  The  analog 
signal  was  passed  through  a  low-pass  filter  whose  cutoff 
frequency  was  less  than  one  half  the  sampling  rate  of 
the  analog  to  digital  converter  (ADC) .  This  procedure 
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assures  that  the  sampling  rate  will  be  at  least  twice 
the  highest  fre  ?ncy  component  of  the  data  and 
minimizes  the  noise  added  by  the  ADC. 

Once  the  data  is  digitized,  it  must  be  converted  to 
units  of  pressure  (psi)  to  be  meaningful.  Program  Slope 
(Appendix  C)  accomplishes  this.  The  program  reads  the 
calibration  information  at  the  beginning  of  the  data, 
asks  the  operator  what  pressure  this  data  represents. 

The  program  then  computes  a  slope  and  intercept,  which 
is  applied  to  each  data  point.  The  particular  data  I 
used  was  inverted  in  the  process  of  digitizing, 
therefore,  this  program  will  invert  the  data  again  to 
return  it  to  its  proper  phase  (Figure  19). 

At  this  point  the  data  must  be  edited  to  remove  all 
signal  prior  to  the  initial  rise  of  the  blast  wave. 

This  editing  is  done  because  program  Filter  assumes  that 
the  blast  wave  begins  at  time  equal  to  zero.  If  the 
editing  is  not  done,  a  phase  shift  will  result  in  the 
filtered  data. 

The  final  step  is  to  actually  filter  the  data  using 
program  Filter  (Appendix  D) .  This  program  asks  the 
operator  for  several  pieces  of  information.  First  it 
asks  for  equivalent  explosive  weight  and  the  distance  to 
ground  zero.  It  must  have  this  information  to  compute 
the  ideal  blast  wave.  The  next  request  the  computer 
makes  is  for  the  power  of  two  for  the  FFT.  The  number 
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input  at  this  point  will  determine  the  resolution  of  the 
frequency  spectrum  computed  for  both  the  ideal  and 
actual  blast  waves.  It  should  be  such  that  two  raised 
to  this  power  is  greater  than  twice  the  number  of  data 
samples  in  the  actual  data  and  at  the  same  time  less 
than  4097  to  prevent  going  outside  the  bounds  of  the 
program. 

The  next  request  of  the  program  is  for  the  operator 
to  input  the  time  per  sample  for  the  actual  data.  This 
is  used  to  correlate  the  ideal  blast  wave  to  the  actual 
data. 

The  final  request  is  to  input  the  desired  percent  of 
peak  overpressure  (95-100).  If  the  response  is  100, 
then  the  Frielander  approximation  is  used,  but  if  95-99 
is  chosen  then  equation  (3)  is  used  with  the  appropriate 
value  selected  for  the  cutoff  frequency.  In  most  cases, 
this  choice  will  make  very  little  difference  in  the 
output  data.  Therefore,  I  recommend  100%  be  used  to 
decrease  the  computation  time  slightly. 

At  this  point,  the  program  computes  the  ideal  blast 
wave  and  its  FFT.  Then  the  FFT  of  the  actual  data  is 
computed.  Since  the  amplitude  of  an  FFT  is  proportional 
to  the  amplitude  of  the  original  signal,  one  can  compare 
the  magnitudes  of  the  two  FFTs.  Comparing  the  magni¬ 
tudes  of  the  FFTs,  the  program  substitutes  the  value  of 
the  FFT  of  the  ideal  blast  wave  for  any  point  in  the  FFT 


that  the  actual  exceeds  the  Ideal  by  2.5%.  This  proce¬ 
dure  has  the  effect  of  removing  only  enough  of  the 
unwanted  frequencies  to  remove  the  noise. 

Program  Filter  requires  that  the  guess  for  the 
explosive  weight  be  higher  than  the  actual  explosive 
weight.  Otherwise,  the  output  data  would  just  be  a 
reproduction  of  the  ideal  blast  wave  for  that  guess.  If 
after  the  first  try  the  data  is  still  noisy,  the  program 
could  be  run  again  with  a  lower  guess  for  the  explosive 
weight. 

The  time  required  for  this  form  of  filtering  is 
going  to  be  dependent  on  the  type  of  computer  used,  but 
should  be  much  faster  than  attempting  to  do  the 
filtering  in  an  analog  fashion.  The  excellent  results 
of  this  method  are  shown  in  Figure  20,  the  filtered 
version  of  Figure  19. 

Although  this  program  reads  the  data  from  a  file  in 
the  computer  and  puts  the  resulting  data  into  another 
file,  there  is  an  underlying  assumption  that  the 
computer  has  the  ability  to  plot  the  data.  Otherwise, 
it  would  be  of  very  little  benefit  to  the  using  agency. 

There  is  always  more  than  one  way  to  accomplish  a 
given  task,  and  so  it  is  in  this  case.  The  following 
section  presents  a  time  domain  method  to  remove  the 
noise. 
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Figure  20.  Data  After  Digital  Filtering 
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5.  Numerical  Smoothing 

2 

A  popular  method  of  smoothing  data  is  the  method 
of  least  squares.  It  is  the  approach  I  have  taken  here. 

Two  conditions  must  be  met  to  compute  an  estimate 
of  data  using  the  method  of  least  squares.  First  there 
cannot  be  any  systematic  errors  in  the  estimate  (i.e.  no 
bias).  For  the  blast  wave  data,  the  average  of  the 
noise  is  zero.  Therefore,  there  is  no  bias. 

The  second  condition  is  that  the  points  computed  as 
a  result  of  the  least  squares  estimation  should  be 
linear  functions  of  the  measurement  data.  This 
condition  creates  a  minor  problem.  The  ideal  blast  wave 
is  not  a  linear  function  of  time.  Therefore,  I  made  a 
slight  change  to  the  standard  least  squares  technique;  I 
call  it  the  "sliding  least  squares." 

Since  the  blast  wave  can  be  approximated  by  a 
series  of  straight  line  segments,  I  wrote  a  program 
which  would  compute  a  sliding  least  square  fit  to  the 
blast  wave  data  (Program  Smooth,  Appendix  E) 

The  program  asks  the  operator  for  the  number  of 
samples  to  "average,”  call  it  (n)  samples.  Then  the 
computer  reads  (n)  samples  from  the  data  file,  computes 
a  standard  least  squares  fit  to  these  points,  and  outputs 
the  first  (n/2)  points  according  to  this  least  square 
estimate. 

2G.  Dahlquist  and  A.  Bjorck,  Numerical  Methods, 
translated  by  N.  Anderson,  (Prentice-Hall,  Inc.,  New 
Jersey,  1974),  pp.  126-129. 
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At  this  point  the  sliding  starts.  The  elements  of 
the  array  holding  the  first  (n)  samples  are  shifted  one 
place  to  the  left.  The  first  data  sample  is  discarded 
and  the  (n  +  1)  sample  is  read  into  the  nth  position  of 
the  data  array.  A  new  least  square  estimate  of  the  data 
is  now  computed  and  used  to  compute  the  (n/2  +  1)  data 
point.  Here  the  process  starts  over  again  by  shifting 
the  data  array  one  place  to  the  left  reading  in  a  new 
point,  etc. 

If  the  end  of  the  data  file  is  reached  before  the 
specified  number  of  points  have  been  calculated,  the 
last  few  points  (as  many  as  n/2)  are  calculated 
according  to  the  latest  least  square  estimate. 

The  result  of  this  method  is  shown  in  Figure  21, 
the  original  test  data  was  that  of  Figure  19. 

Once  the  noise  is  removed,  standard  methods  can  be 
used  to  determine  the  positive  phase  integral  and  the 
peak  amplitude.  The  actual  explosive  equivalence  of  the 
charge  being  tested  can  be  determined  from  these  two 
parameters. 

If  the  data  is  digitized  as  it  is  being  collected, 
the  process  of  smoothing  could  be  done  in  about  the  same 
amount  of  time  as  it  takes  to  convert  the  data  to 
engineering  units.  This  means  a  digital  data 
acquisition  system. 


DIGITAL  TELEMETRY  SYSTEM 


An  analog  system  has  several  disadvantages.  First 
is  the  limit  on  the  number  of  channels  of  data  that  can 
be  recorded  (the  number  of  channels  on  a  tape  recorder) . 
Second  is  the  fact  that  two  power  wires  and  a  coax  cable 
must  be  run  to  each  transducer,  which  means  that  once 
the  transducers  are  located,  they  are  pretty  well 
permanent.  Third,  the  data  must  be  digitized  before  it 
can  be  analyzed  by  a  computer.  The  proposed  system  will 
eliminate  these  restrictions. 

It  is  a  system  using  one  Central  Control  Unit  (CCU) 
and  many  Field  Units  (FU) .  The  number  of  channels  will 
only  be  limited  by  the  number  of  transducers  that  each 
FU  can  service  and  the  size  of  computer  word  used  to 
address  the  individual  FUs.  In  general  the  CCU  will  be 
a  16  bit  computer  capable  of  addressing  up  to  1000  FUs 
via  radio  signals  (2-3  GHz).  Each  FU  must  be  able  to 
service  at  least  eight  (8)  transducers  with  an  overall 
throughput  rate  of  50  000  data  samples  per  second. 

1.  Central  Control  Unit 

This  unit  acts  as  the  "brain”  of  the  system.  Prior 
to  the  test  it  will  automatically  check  all  of  the  FUs 


to  see  if  they  and  their  peripherals  are  operational  and 
ready.  When  the  tine  cones  for  the  test,  the  CCD  will 
tell  the  FUs  to  start  recording  data  at  the  proper 
tines.  Further,  the  CCD  will  start  the  caneras  and 
even  detonate  the  test  item  through  the  FDs.  After  the 
test  is  over,  the  CCD  will  instruct  each  of  the  FDs,  one 
at  a  tine,  to  transmit  the  recorded  data  back  to  the 
main  computer  for  storage  and  analysis  at  a  later  tine. 
Finally,  the  CCD  will  make  a  post  operational  check  of 
each  of  the  FDs  and  "flag”  any  unit  that  is  question¬ 
able,  thereby,  allowing  that  unit  to  be  checked  immedi- 
ately  for  the  cause  of  the  problem  and  a  possible  data 
bias.  At  this  point,  the  data  could  either  be  recorded 
on  tape  and  sent  to  the  main  computer  on  base  for  anal¬ 
ysis,  or  it  could  be  digitally  filtered  as  discused 
above  and  then  recorded  and  sent  to  the  base  for  further 
analysis. 

The  CCD  will  consist  of  an  electronic,  digital 
computer  (minicomputer  or  main  frame)  and  a  radio 
transceiver  operating  in  the  range  of  2-3  GHz  connected 
by  the  appropriate  interface. 

The  computer  shall  be  programable  in  a  high  level 
language  (Fortran,  Algol,  or  Pascal),  with  instructions 
which  permit  the  addressing  of  peripherals  from  within 
the  program.  As  a  minimum  this  computer  must  have  the 
computing  power  and  internal  memory  (one  megabyte) 
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equivalent  to  an  HP3000/33,  with  some  form  of  graphics 
capability  (i.e.  be  able  to  make  two-dimensional ,  hard 
copy  plots  of  data) . 

The  following  peripheral  equipment  is  required: 

1.  One  (1)  multisurface  disc  memory,  with  expansion 
capability  to  five  (5)  such  units.  Each  unit  shall  have 
a  minimum  of  120  megabytes  of  disk  storage. 

2.  Industry  standard,  IBM  compatable,  magnetic  tape 
drive,  nine  (9)  track,  capable  of  reading  and  writing 
1600  bits  per  inch. 

3.  A  video  type  terminal  with  graphics. 

4.  Lineprinter. 

5.  Graphics  printer. 

6.  Transceiver. 

7.  Card  reader. 

The  above  items  could  be  combined,  for  example; 
items  3,  4,  and  5  could  be  provided  in  one  unit,  the  HP 
2648A  Graphics  Terminal/Printer. 

In  addition  to  the  above,  a  minimum  of  five  (5) 
extra  peripheral  ports  must  be  included  with  at  least 
two  (2)  being  serial  RS-232C  compatible  and  at  least  one 
(1)  shall  be  a  16  bit  parallel  port.  These  would  be 
available  for  future  expansion  of  the  system. 

The  computer  shall  have  a  diagnostic  program  to  aid 
in  trouble-shooting  repair  problems.  This  program  will 
diagnose  problems  down  to  the  "board”  level,  internally. 


and  indicate  the  peripheral  port  in  the  case  of  external 
problems. 

One  further  requirement  of  the  CCU  is  that  it  be 
able  to  address,  through  the  transceiver,  each  of  the  FOs 
individually  as  well  as  collectively,  individually  to 
allow  for  checking  of  status  and  to  transmit  data,  and 
collectively  for  starting  time  messages  etc. 

The  transceiver  interface  to  the  computer  shall  be 
designed  to  accept  transceivers  manufactured  by  at  least 
two  (2)  different  companies  and  the  interface  shall  be 
independent  of  the  type  transmission  used  by  the 
transceiver. 

2.  Transceivers 

The  exact  format  for  the  transceiver  shall  be  left 
to  the  discretion  of  the  contractor.  It  will  be  a  form 
of  Frequency  Modulation  with  some  arrangement  to  insure 
positive  communication  is  established  before  significant 
information  is  transmitted.  The  framing  arrangement  for 
the  data  transmission  will  likewise  be  left  to  the 
discretion  of  the  contractor.  The  minimum  acceptable 
range  for  the  telemetry  system  is  25  miles. 

Application  will  be  made  to  the  Federal  Communica¬ 
tion  Commission  for  a  specific  frequency,  bandwidth,  and 
power  allocation  just  prior  to  releasing  the  contract 
for  bidB. 


The  transceiver  in  the  PD  will  obviously  have  to  be 
compatible  with  the  one  in  the  CCD.  The  FD  transceiver 
must  occupy  no  more  than  216  cubic  inches  and  be  battery 
powered.  The  CCD  will  use  the  most  efficient  antenna 
possible ,  while  each  of  the  FDs  will  use  only  the 
minimum  possible  antenna  to  limit  the  area  they  present 
to  the  blast  wave  and  also  limit  the  cost.  As  presently 
planned,  the  CCD  will  be  located  in  line-of-sight  with 
the  FDs  and  slightly  higher  than  the  FDs.  If  the  test 
site  location  is  ever  changed,  a  repeater  will  probably 
have  to  be  installed  to  handle  the  transmission. 

3.  Field  Onit 

Each  FD  is  a  totally  independent  unit  and  operates 
entirely  on  its  own  once  a  start  signal  is  received. 

The  FU  must  be  able  to  check  itself  and  the  transducers 
connected  to  it,  and  then  transmit  a  status  report  of 
both  to  the  CCD.  It  must  be  capable  of  determining  when 
to  start  and  stop  recording  data.  The  FD  must  be 
capable  of  sampling  and  recording  data  at  rates  up  to 
and  including  50  000  samples  per  second,  overall 
throughput,  from  programably  selectable  channels. 

When  the  signal  is  received  from  the  CCD  that  the 
test  item  will  detonate  in  a  given  amount  of  time, 
several  events  must  take  place  in  the  FD.  First  a  real¬ 
time  countdown  clock  is  started.  Second,  two  calibration 
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values  are  read  for  each  transducer  and  recorded  as  the 
first  two  pieces  of  data  for  each  transducer.  These 
calibration  values  are  to  be  a  voltage  substition  type, 
with  the  calibration  voltage  switched  into  the  system 
prior  to  the  first  amplification.  One  voltage  repre¬ 
sents  the  highest  expected  output  of  the  transducer  the 
other  is  to  be  zero  volts.  These  calibration  values 
must  be  accurate  to  within  0.1%  of  the  desired  value. 
Third,  a  relay  closure  (1.5  amperes  max,  120  volts  ac) 
must  occur  at  a  preprogrammed  time,  plus  or  minus  0.5 
milliseconds. 

The  ability  of  the  FU  to  provide  a  relay  closure, 
at  a  precise  time  prior  to  detonation,  is  for  starting 
cameras  or  even  detonation  of  the  test  item. 

The  fourth  event  will  be  to  start  recording  data 
from  the  preprogramed  number  of  transducers  at  the 
appropriate  sampling  rate.  This  data  will  be  rotated  in 
memory  until  actual  data  (rise  of  the  blast  wave)  is 
detected.  Then  the  data  10  milliseconds  prior  to  the 
rise  and  240  milliseconds  after  the  rise  will  be 
recorded  and  stored.  Also  stored  for  each  transducer 
will  be  the  real-time  clock  word  corresponding  to  the 
first  data  sample. 

The  fifth  event  is  a  post-calibration,  to  be 
performed  identical  to  the  first  one.  With  these  two 
pieces  of  information  stored  as  the  last  two  words  in 


the  data  file.  This  post-calibration  should  take  place 
at  least  30  seconds  after  the  last  data  sample  is 
acquired. 

The  final  event  the  FU  must  accomplish  on  its  own 
control  is  to  output  the  data  to  a  cassette  recorder. 
This  could  be  started  as  soon  as  the  post-calibration  is 
complete.  The  data  is  to  be  put  on  the  cassette  one 
transducer  at  a  time  with  enough  leader  information  to 
identify  the  transducer  location,  time,  and  date  of  the 
test. 

The  last  event  of  the  test  would  be  for  the  CCU  to 
sequentially  have  each  of  the  FOs  transmit  the  data  they 
collected  back  to  the  CCD  in  the  same  format  as  the  FU 
used  to  record  the  data  on  cassettes. 

The  FU  should  be  versatile  enough  that  any  of 
several  types  of  transducers  could  be  connected  with 
little  or  no  modification  to  the  FU.  The  FU  must  be 
battery  operated  to  eliminate  the  need  to  run  any  power 
lines  or  signal  lines  at  all  except  to  the  transducers 
it  services.  This  battery  must  be  capable  of  operating 
the  FU  continuously  for  up  to  four  (4)  hours. 

The  FU  consists  of  a  transceiver,  connected  via  an 
appropriate  interface  to  the  microprocessor.  The 
microprcessor  controls  all  the  functions  mentioned  above 
and  acts  as  a  traffic  controller  to  pass  data  from 
either  method  of  storage  to  the  transceiver  for 
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transmission  to  the  CCU.  The  microprocessor  will  be  16 
bit  with  a  sufficient  instruction  set  to  accomplish  the 
above  described  tasks. 

The  random  access  memory  (RAM)  shall  be  sufficient 
to  store  12  500,  16  bit,  samples  plus  the  necessary 
program  Instructions.  If  the  sampling  rate  is  50  000 
samples  per  second,  this  represents  250  milliseconds  of 
signal  data,  which  is  twice  that  of  most  of  the  signals 
measured  during  the  testing  of  1980.  The  erasable, 
programable,  read-only  memory  (EPROM)  is  used  to  store 
the  program  so  it  can  easily  be  changed  to  accommodate 
changes  in  use  of  the  FU.  For  example,  suppose  digital 
signals  were  being  received  from  several  different 
velocity  transducers  instead  of  the  digitized  signal 
from  the  ADC,  a  different  format  might  be  required  to 
map  the  data  into  memory.  Definitely,  there  would  be  big 
differences  in  the  self-test  function  and  therefore,  a 
modified  program  would  be  required. 

Of  the  channels  each  FU  sevices,  at  least  four  (4) 
will  be  analog  and  four  (4)  will  be  digital.  The  digital 
channels  will  be  16  bit  parallel  ports. 

Each  of  the  analog  channels  will  have  an  adjustable 
gain  amplifier  (1-1000)  with  a  selectable  input  impedance 
of  50  ohms  or  10  megohms,  also  to  be  selectable  is  single 
or  differential  input.  Each  of  the  analog  channels  will 
incorporate  a  single-pole,  low-pass  filter  whose  cutoff 


frequency  (-3  dB)  is  continuously  adjustable  from  100  Hz 
to  20  kHz.  The  amplifier  must  have  a  constant  gain  down 
to  and  including  dc. 

The  multiplexer  must  be  able  to  select  any 
combination  of  the  analog  channels  with  uniform  sampling 
over  those  channels  selected. 

The  ADC  must  accept  a  full  scale  input  of  plus  to 
minus  10  volts  and  provide  at  least  12  bit  resolution, 
including  sign,  at  the  output.  This  is  to  be  accom¬ 
plished  at  a  minimum  throughput  rate  of  50  000  samples 
per  second. 

The  data  acquisition  section  of  the  FUs  will  meet  or 
exceed  the  following  specifications: 

Accuracy:  0.025% 

Stability:  0.025% 

Temperature  Coefficient:  0.01%/°C  (from  0-70  °C) 
Linearity:  0.025% 

Repeatablity:  0.025% 

Sample  and  Hold:  <  100  nanosecond 

aperature  uncertainty 

Output  code:  Binary 
Crosstalk:  <  0.025% 

Common  Mode  Rejection:  >75  dB 

Noise,  3  sigma:  <600  microvolts,  peak-to-peak 
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The  interface  between  the  microprocessor  and  the 
transceiver  must  be  similar  to  the  one  specified  for  the 
CCU  in  that  it  must  accept  transceivers  manufactured  by 
at  least  two  different  companies. 

The  EMP  of  the  explosion  can  be  quite  large, 
depending  on  the  Bize  of  the  explosion  and  distance  from 
ground  zero.  Therefore,  an  additional  requirement  of 
the  FU  will  be  to  electically  disconnect  the  antenna 
from  the  front-end  of  the  transceiver  one  second  prior 
to  detonation  and  reconnect  it  about  one  minute  after 
the  detonation.  This  will  help  prevent  the  transceiver 
from  being  damaged  by  the  EMP.  Both  the  antenna  and  the 
transceiver  front-end  should  be  connected  to  ground 
during  this  period. 

The  final  requirement  of  the  FU  is  that  it  fit  in  a 
box  2  ft  x  2  ft  x  2  ft.  This  box  will  be  buried  in  the 
ground  with  the  top  of  the  box  flush  with  the  ground. 

This  box,  therefore,  will  have  to  be  air  and  water 
tight.  The  antenna  can  be  mounted  on  the  top  of  the  box 
and  feed-throughs  for  the  transducers  can  be  mounted  on 
the  sides  of  the  box.  The  FU  must  be  able  to  operate  in 
this  environment  and  it  must  be  easily  removeable  from 
the  box  (i.e.  disconnect  antenna  and  transducers,  lift 
the  FU  from  the  box). 

The  initial  cost  of  this  system  may  be  more  than  a 
comparable  analog  system,  however,  the  versitility  and 
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portability  make  the  telemetry  system  much  more  use- 
able.  Further,  in  the  long  run  this  system  will  save 
money,  because  the  cost  to  change  transducer  locations 
should  be  cut  by  a  factor  of  ten,  other  savings  will  be 
realized  in  data  reduction. 


l 


RECOMMENDATIONS 


To  produce  data  that  will  be  readily  usable  by  the 
requesting  agency,  I  recommend  that  the  above  specified 
digital  acquisition  system  be  purchased  and  used  when 
testing  explosive  items.  Further,  the  use  of  Figure  6 
and  Appendix  F  to  set  the  cutoff  frequency  of  the  low- 
pass  filter  will  help  reduce  the  noise  recorded  and  also 
allow  the  actual  peak  overpressures  to  be  calculated 
from  the  recorded  data. 

Once  the  data  is  recorded,  it  should  be  filtered 
with  one  of  the  two  programs  provided  in  Appendices  D 
and  E.  This  will  give  the  using  agency  the  most  ideal 
data  possible  to  use  in  their  calculations  and  decisions. 


APPENDIX  A 


PROGRAM  RESPONSE 
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progran  RESPONSE 

c*  This  progran  was  developed  as  part  of  an  Electrical 
c*  Engineering  Thesis  by  Theodore  J  Moody,  entitled, 
c*  "Instrunenting  an  Explosive  Test  Arena,"  at  the 
c*  University  of  Utah,  May  1982. 

c ******************************************************* 

♦control  nolist 
♦control  f ile=25-39 
i2p  =  9.0 
n  =  2*#i2p 

c*  Operator  inputs  a  guess  for  the  pole  location 

display  'input  pole  Multiplier' 
accept  pnult 
1  =  0 

do  1500  11  =  1,4 

display  '  Lanbda  LPF  cutoff  Hz  Peak' 
do  1000  12  =  1,50 
1=1  +  1 

c************************************ 

c*  Lanbda  =  range/ (explosi ve  weight  A(i/3)> 

c *  where  xlan  =  lanbda 

c*  r  =  range 

c#  wt  =  explosive  weight 

c*  This  paraneter  is  fron  The  Air  Force  Technical 

c*  ilA-i-47,  Departnent  of  Defense  Explosive  Hazard 

c#  Classification  Procedures,  March  1981 . 

c ************************************ 

xlan  =  10 . 0**( <0 . 0165*(f loat ( 1) ) )-0 .301 ) 
c ************************************ 

c*  The  peak  overpressure  and  inpulse  described  in 

c*  in  this  TO  are  only  described  for  0 . 5< lanbda< 1000 . 

c *********************** ************* 
c ************************************ 

c*  The  inpulse  of  the  peak  overpressure  is  known  only 
c*  exper inen tially ,  therefore,  two  equations  are 

c*  needed  to  decribe  the  positive  inpulse;  one  for 

c*  the  range  0 . 5<lanbda<2 . o48  (the  equation  following 

c*  label  410)  and  one  for  the  range  2. 048< lanbda< 1000 

c*  (the  equation  with  label  420). 

c*  xinp  *  Positive  inpulse 

c#  The  following  calculates  the  value  of  the  inpulse 
ct  for  the  given  lanbda . 
c #**##**####*******#**#****#****##*#* 
xi  =  40.0 
y 1  =  0.05 
x2  =  2.048 
y2  *  11.474 
x3  =  1.442 
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y3  *  13.16 

x4  =  1.0 

y4  =  20.0 

xll  =  alogl0<xi> 

yll  =  alogiO(yl) 

x21  =  aloglO<  x2> 

y21  =  aloglO(y2) 

x31  =  aloglO(x3> 

y31  =  aloglO(y3) 

x41  =  alogl0<x4> 

y41  =  alog!0(y4) 

a  *  <yil-y21)/<xll-x21> 

b  =  yll  -  a*xil 

b  =  10.0**b 

c  =  ( <  y31-y41 )*<  x21-x41 >-<  x31~x41 )#<  y21-y41 >  >/ 

1  ( <  <x31*x31-x41*x41>*<x21-x41>>-<  <x31-x41>* 

2  <x21*x21-x41*x41> > ) 

d  =  <y21-y41-c*<x21*x21-x41*x41>)/<x21-x41> 
e  *  <y41-<c*x41*x41>-d*x41> 

410  if ( xlan . ge . 2 . 048)go  to  420 
c******************************************************** 


c*  Calculates  impulse  if  in  parabolic  portion 

f  =  aloglO(xlan) 
g  *  <c*f*f >+<d*f >+e 
xinp  =  <10.0**g) 

c*  Calculates  impulse  if  in  linear  portion 
go  to  430 

420  xinp  =  b*( ( xlan)**a ) 

430  continue 

c************************************ 


c*  Sends  lambda  and  weight  to  files  number  31  and  32 

c************************************ 
wt  =  (1 . 0/xlan>**3. 0 
wr ite(31 ) xlan 
write(32>xinp 

c ******************************************************** 
c*  The  following  equations  are  from  an  article  by 

c*  L.  Ciglio-Tos  and  T.E.  Linnenbrink,  “Air  Blast 

c*  Pressure  Measurement  Systems  and  Techniques", 
c*  MINUTES,  15th  DODESB  Seminar,  September  1973, 
c*  p.  1400.  Although  the  equations  in  the  reference 
c#  have  considerable  overlap,  Lambda  »  71.725155  was 

c*  choosen  as  the  change  over  point  because  this  is 

c*  the  point  where  the  two  equations  are  equal  (to 
c*  9  decimals > . 

c *  The  first  equation  calculates  the  peak  overpressure 

c*  <pm>  in  the  range  0 . 5< pm<71 . 725155 

c ******************************************************** 

if (xl an. gt. 71 . 725155) go  to  440 


xlnpn  =  7.0452041-1 . 6277561* ( al og (  x lan) )  -  0 . 27399088* 

1  ( (alog< xlan) >**2) 

2  -0 . 065973136*< (alog < xlan) >**3>+0 . 0065412563* 

3  < <alog<xlan) >**4> 

4  +0 . 04B236359*< (alog(xlaM) >**5)-0 . 020072553* 

5  ( <al og < xlan ) >**6> 

6  +0  003019Q449*<<alog<xlan>>**7>-0. 00015984026* 

7  ( (alog ( xlan) )**8) 
pn  *  exp(xlnpn) 

go  to  450 

^  |||  ||^  ^  ^  ^  |||  j||  j||  ^  ^  ^  ^  ^ 

c*  Calculates  peak  overpressure  <pn)  in  the  range 
c*  71 . 725155<pn< 1000.0 

440  p«  =  226. 61762/< <xlaM)**(i . 4065913)) 

c*  Calculates  the  tine  of  the  positive  phase  of  the 

c*  blast  wave  <xto>  froM  an  analytical  integral 

c*  solution  of  the  ideal  blast  wave. 

o«  i  i  i  i  3k  )k  *  Jk  3k  j|[  sk  3k  Jk  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k  3k 
450  x t o  ■  0 . 001*xiMp*2. 718281828/pM 
c ************************************ 

c*  This  equation  is  derived  fro«  solving  the  integral 

c*  of  p<t>  =  pn*exp (~t/t o >*( i-t/t o ) 

c*  the  Frielander  Appr oxination  to  a  Blast  Uave,  in 

c*  combination  with  the  above  equations  for  pn.  The 

c*  object  is  to  find  the  tine  the  overpressure 

c*  returns  to  anbient  (xto). 

c ************************************ 
write(33>xto 
wr ite(34)pn 

c*  The  following  uses  equation  (3>  developed  in  the  Main 
c*  part  of  the  theses  Section  A 

c*  This  calculates  the  value  of  y(t)  until  a  Maxinun  is 
c#  detected,  then  outputs  the  naxinun  and  calculates  the 
c*  Range-Bandwid th  Product 

to=xto 

pi  -  3.141593 
ainc  *  pi/to 
a  *0.0 

tine  *  to/10000.0 

a  *  pnult*ainc 

ato  *  a#to/<a*to-i . 0) 

t  «  0.0 

ynnl  *  -1.0 

do  100  i  *  1,9000 

t  *  t+tinc 

b  *  -t/to 

c  *  exp(-a*t) 

d  *  exp(b) 
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Program  SPECTRUM 

c *********************************************** ******** 

c*  This  program  was  developed  as  part  of  an  Electrical 
c*  Engineering  Thesis  by  Theodore  J  Moody,  entitled, 
c*  "Instrumenting  an  Explosive  Test  Arena,"  at  the 
c*  University  of  Utah,  May  1982. 

c ******************************************************* 
c ************************************ 

c*  This  program  reads  a  file,  takes  the  Fourier  Transform 
c*  of  that  file  and  outputs  the  Magnitude  of  the  transform 
c*  and  the  Power  Spectrum, 
c ******* ********* ******************** 

♦control  nolist 
♦control  file=25-29 

dimension  u(4096),v(4096),w(4096) 
c ********************************* 

c*  Reads  data  from  file  ftn25  designated  by  the  operator 
c********************************* 
do  100  k  *  1,4096 
read<25,end=120)u (k> 
v ( k )  =  0.0 
100  continue 

display  '  Read  data  ' 

120  k  =  k-i 

c*  Operator  selects  the  power  of  two  for  FFT 
c ************************************ 

125  wr ite<6 , 130 ) 

130  f ormat < lx ,24hPower  of  2  for  the  FFT  => 

read  <5,*) i2p 
140  format ( 15) 

n  =  2**i2p 
nm  =n-l 

if<n-k >142,155,145 

142  write<6,143> 

143  format < lx ,20hPower  of  2  too  small) 

go  to  125 

145  if <n . It ,4097)go  to  147 

write<6, 146> 

146  format < 1 x , 20hP ower  of  2  too  large) 

go  to  125 

147  do  150  j  «  <k+l),n 
u(j)  =0.0 

v ( j>  =0.0 
150  continue 

155  continue 

pi  *  3 . 141592654 

c******************************** 

c*  Developes  window  to  properly  condition  the  data 
c*  before  transforming.  Uses  the  Blackman  window 
c******************************** 
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do  170  j  =  i,n 

w<  3 )  =  0 . 42- <  0 .5*cos<2*pi*< j-l>/nM>> 
i  +(0 ,08*cos<4*pi*< j~l)/nm>> 

170  continue 

display  '  Window  built  ' 
ns  =  n/2 

C**#*#**##*##******##*JM***#**#* 

c*  Shifts  window  over  data 
c  *  ##**#**#**#  *  *#  #*  *#*#  *  ik  *#  *#**  *  * 

call  shif t (w, n ,ns) 

c*  Windows  data 

c  *##*####  ##  *##  *  *#  *  **  %  *#*  *  *  *#*  *** * 

do  180  j  *1,0 
uCj)  =  w<  j )#u<  j  > 

180  continue 

display  '  data  windowed' 
c* fc***#*#*#**#*****##*#*****#* 
c#  Takes  Fourier  Transform 
c***************************** 
call  nuf f t < u , v , i2p , 1 > 
do  185  j=i,n 
wr  ite(28 ) u  <  j  > 
write<29)v<j) 

185  continue 

display  '  FFT  filed' 
endfile  28 
endfile  29 
ymax  *  0.0 
ymin  *  1000.0 
znax  -0.0 
zmin  *  10000 

c * ** * * * # * * M * * * * * * * * * * * * * * * ##* * # * * 
c#  Determines  Magnitude  of  transform 

C  ******* ************************* 

do  190  j  *  1,  n 

u(j>  *  sqr t < < u< j )**2 )+ < v < j )#*2) ) 

c ******************************* 

c*  Determines  the  Maximum  and  minimum  values  of 
[  c#  the  magnitude 

c ## * * * #* * * * * * * * * * ** * * # * * * * * * * * # * 

if <u< j) ,gt .yMax)ymax  *  u(j) 
if <u<3> . It ,ymin>ymin  *  u<j> 
c ****************************** 
c*  Outputs  Magnitude  to  file  ftn26 
c ****************************** 
write<26>u<  j> 

c********************************* 

c#  Determines  the  Power  Spectrum  of  the  transform 
c********************************* 
v(j)  *10 . 0#a log  10 ( u < j ) > 
c *********************************** 
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cHe  Determines  the  naxinun  and  Minimm  values  of 
c*  the  Power  SpectruM 
c He#**# ************ * He* ********** **#* * 
if  <v< j) .gt .z«ax)zMax  =  v(j) 
if  ( v ( 3 ) . It . znin) zhin  =  v<j) 

^  ijc  i  i  He  He  ijt  i  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He 
cHe  Outputs  the  Power  SpectruM  to  file  ftn27 
cH(  He  He  He  He  He  He  He  He  He  He  He  H<  He  He  He  He  He  He  He  He  He  He  He  He  Hi  He  He  He  He  He  He  He  He  He  He  He  He  He 
wr i te< 27) v  <  j ) 

190  continue 


display  '  Magnitude  and  Power  SpectruM  Filed' 


cHe  Outputs  MaxiMun  and  MiniMUM  values  of  the  Magnitude 
c#  and  the  Power  SpectruM  to  the  teletype 


200 
2i  0 
22G 
230 


write(fe.200)yrtax 
wr ite<6,210 ) y«in 
write(6,220)z«ax 
write(6,230)zMin 

f oraat ( lx ,24hFreq  Magnitude  MAXIMUM 
f ornat ( ix >24hFreq  Magnitude  MINIMUM 
f ornat ( ix >24hPower  SpectruM  MAXIMUN 
f ornat < ix ,24hPower  SpectruM  MINIMUM 
endfile  26 


= ,f 12 . 5) 
- ,f 12 . S) 
= ,f 12 . 5) 
= ,  f  12  .  S ) 


endfile  27 


stop 

end 


CHe 

This  is  a 

subroutine  to  perf< 

irn 

the 

Fast 

CHe 

Fourier  Transform  (FFT)  of  a 

c  onplex 

array 

CHe 

where : 

CHe 

XR 

=  Real  array 

C* 

XI 

=  Imaginary  Array 

CHe 

M 

=  The  order  of  the 

FFT 

<2A 

M) 

CHe 

IDIRN  =  The  direction 

of 

the 

Tr ansf orn 

CHe 

1  =  forward 

CHe 

-1  =  reverse 

CHe 

The  result 

returns  in  XR  and 

XI 

C  *  He  H<  He  He  He  He  Hi  He  He  He  He  He  He  He  He  He  Hi  He  He  He  Hi  He  Hi  He  He  He  He  He  He  He  He  He  He  He  He  He  He 
SUBROUTINE  NUFFT C XR , XI ,M , IDIRN > 
DIMENSION  XR ( 1 ) , XI ( 1 > 

N=2**M 
PI*3 . 141S9H7 
DO  20  L=1,M 
LE=2*He<M+i-L> 

LEl=LE/2 
UR  =  i  .0 
UI«0.0 

WR=COS( PI/FLOAT (LEI ) ) 

WI=-SIN< PI /FLOAT <LE1) ) 

DO  20  J=i,LEi 
DO  10  I=J,N,LE 
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IP=I+LEi 
TR=XR ( I >+XR < IP  > 

TI=XI<I)+XI<IP> 

TMR=XR ( I ) -XR ( IP  > 

TMI=XI(I)-XI<IP> 

XR  < IP)=TMR*UR-TMI*UI 
XI < IP  >=TMR*UI+TMI*UR 
XR  < I >=TR 
10  XI < I ) “TI 

TR=UR#WR-UI*WI 
UI=UR*WI+UI*WR 
20  UR=TR 

NV2=N/2 

NMi-N-i 

J=i 

DO  30  1=1, NMi 

IF ( I  .GE.  J)  GOTO  25 

TR=XR< J) 

TI=XI <  J) 

XR  <  J ) =XR ( I ) 

XI(J)=XI<I> 

XR ( I )=TR 
XI(I)=TI 

25  K=NV2 

26  IF < K  GE.  J)  GOTO  30 

J=J-K 

K=K/2 
GOTO  26 
30  J=J+K 

IFdDIRN  .EQ.  1>  GOTO  77 
DO  99  1=1, N/2 
TTR=XR(I)/N 
TTI=XI(I)/N 
XR  < I )=XR (N-I  +  l )/N 
XI ( I >  =XI ( N-I  +  l > /N 
XR (N-I+l )=TTR 
XI (N-I+i )=TTI 
99  CONTINUE 

TTTR  =  XR  <  N ) 

TTTI=XI (N> 

DO  777  I=N,2 ,-i 
XR ( I )=XR  < I — 1 ) 

XI<I>-XI<I-i> 

777  CONTINUE 

XR  < 1 >=TTTR 
XI  < 1 ) “TTTI 
77  RETURN 

END 

CftM**************************************************** 
C*  THIS  SUBROUTINE  SHIFTS  THE  ARRAY  TO  THE  RIGHT 

C*  N/2  SAMPLES  IN  A  CIRCULAR  SHIFT 

C******************************************************* 


SUBROUTINE  SHIFT < X ,LS ,NS> 
DIMENSION  X(LS> 

DO  200  J  *  1  ,NS 
Y  ■  X(i) 

DO  100  K  =  2 , LS 
X(K-i)  *  X  <  K  ) 

X ( LS )  »  Y 
RETURN 


APPENDIX  C 
PROGRAM  SLOPE 


1 
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Program  SLOPE 

c******************************************************* 
c*  This  program  was  developed  as  part  of  an  Electrical 
c*  Engineering  Thesis  by  Theodore  J  Moody)  entitled, 
c*  "Instr umenting  an  Explosive  Test  Arena,"  at  the 
c*  University  of  Utah,  May  1982. 

c******************************************************* 
c *  This  program  reads  data  fro*  a  file  and  converts 

c#  the  data  to  a  units  of  pressure.  It  reads  the 

c#  calibration  information  at  the  beginning  of  the 

c*  data,  asks  the  operator  what  pressure  this  repre- 

c#  sents,  then  computes  a  slope  and  intercept,  which 

c*  is  applied  to  each  data  point.  Also  inverts  the 

c*  data. 

♦control  nolist 
♦control  file  =  25-35 

integer  rfln,wfln 

display  'input  file  number  to  work  on' 
accept  ifln 

display  'number  of  points  in  top  step' 
accept  nt 

display  'starting  point  for  top  step' 
accept  nts 

display  'number  of  points  in  zero' 
accept  nz 

display  'starting  point  for  ZERO' 
accept  nzs 

c ******************************************************** 
c*  Reads  TOP  STEP  information  and  averages  it 

do  100  j  =  i , <nts~i ) 

100  read < if In ,*) x 
sum  *0.0 
k  =  0 

do  200  j  *  nts, <nts+nt-i) 
read(ifln,*,end=300)x 
sum  *  sum  -  x 
k  *  k  +  1 
200  continue 
300  x2  =  sum/k 
rewind  ifln 
sum  *0.0 
k  *  0 

c ******************************************************** 
c*  Reads  ZERO  information  and  averages  it 

do  400  j  e  i,<nzs-l> 

400  read( if In , *> x 

do  500  j  *  nzs, <nzs+nz-l > 
read (if In ,*,end*600 ) x 
sum  *  sum  -  x 


u  u 


73 


k  ■  k  +  i 
500  continue 
600  xi  »  sum  /  k 

display  'TOP-STEP  represents  what  pressure' 
accept  y2 

c*  CoMputes  Slope  and  Intercept  (xicp)  of  the  data 

^  ^  j||  ^  ^  ^  ^  ^  j|^  ^  ^  ^  j||  ^  ^  ^  ^  ^  j|^  ^  ^  ^  ^  ^  ^  ||^  ^  ^  ^  ^  j|^  ^  ||^  ^  ^ 


slope  as  y2/(x2-xl> 
xicp  as  -xi  *  slope 
display  'Slope  =  ', slope 
display  'Intercept  =',xicp 
30  f orMat ( ix , 16, 5x , gi2 . 5> 

•40  f  orMat  (  lx ,gi2 . 5) 

C**************************** !*#********#****#*********#** 

c*  Applies  slope  and  intercept  to  file  to  produce  a 
c*  data  file  of  the  desired  length  in  engineering  units 

display  'how  Many  points  in  current  file' 
accept  i 

50  display  'how  Many  points  in  finished  file' 
accept  n 

if (  i  .  le  .  n  )go  to  iiO 

display  'finished  file  needs  to  be  BIGGER!!!!!' 
go  to  50 
iiO  k  =  0 


of  file  to  read  fron' 


display  'input  slope' 
accept  slope 

c  display  'input  intercept' 

c  accept  xicp 

display  'Input  file  nuMber 
accept  rfln 

display  'Input  file  nuMber  of  file  to  write  to' 
accept  wfln 

c ******************************************************** 

c*  Applies  slope  and  intercept 

c ******************************************************** 
do  210  j  =  l,i 
readtrfln ,*,end=310>x 
k  *  k  +  1 


x  *  (-x)*slope  +  xicp 
wri te< wf In , 40  >  x 
210  write<6,30)k,x 
310  if(k.ge.n)go  to  1000 
do  410  j  ®  <k+i),n 
410  write(wfln,40) 0 . 0 
display 

1  <n-k-l ),' zeros  were  added  to  file', wfln 
go  to  1100 

1000  display  'There  were  no  zeros  added  to  the  file' 
1100  stop 
end 


APPENDIX  D 
PROGRAM  FILTER 


Program  FILTER 

c******************************************************* 
c*  This  program  was  developed  as  part  of  an  Electrical 
c#  Engineering  Thesis  by  Theodore  J  Moody,  entitled, 
c*  “Instrumenting  an  Explosive  Test  Arena,"  at  the 
c#  University  of  Utah,  May  1982. 

c******************************************************* 
♦control  nolist 
♦control  file  =  25-35 

dimension  x <4096 >, y (4096) , z (4096) 
integer  f, percent 
pii  =  3  141593 

c ******************************************************* 
c*  Operator  inputs  guess  for  explosive  weight,  should 
c*  be  larger  than  actual 

c ******************************************************* 

wr i te<  6 , 10  > 

10  f orma t < 1 x , 33h 

1  Input  equivalent  explosive  weight) 

accept  wt 
wr ite<6 ,20 ) 

20  f ormat ( lx , 29h 

1  Input  distance  to  ground  zero) 

accept  range 
r  *=  range 

c************************************ 
c*  Operator  selects  the  power  of  two  for  FFT 

125  write<6,130> 

130  f ormat < lx ,24h Power  of  2  for  the  FFT  => 

read  <5 , #) i2p 

140  formate  IS) 

n  =  2**i2p 

if <n .It  4097)go  to  147 
wr ite  <  6 , 146) 

146  format < ix ,20hPower  of  2  too  large) 

go  to  125 

147  xlam  *  r/ < wt**< 1 . 0/3 . 0 ) ) 

c ************************************ 

c*  The  impulse  of  the  peak  overpressure  is  known  only 
c*  experimentally ,  therefore,  two  equations  are 
c*  needed  to  decribe  the  positive  impulse)  one  for 

c*  the  range  0 .5< lambda<2 . o48  (the  equation  following 

c*  label  410)  and  one  for  the  range  2 . 048< lambda< 1000 
c*  (the  equation  with  label  420). 

c*  ximp  *  Positive  impulse 

c*  The  following  calculates  the  value  of  the  impulse 

c*  for  the  given  lambda, 

c************************************ 
xi  *  40.0 
yl  **  0.05 
x2  «  2.048 
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y2  =  11,474 

x3  -  1.442 

y3  *  13.16 

x4  *  1.0 

y4  =  20.0 

x  1 1  =  aloglO <  xl ) 

yll  =  aloglO(yl) 

x21  =  alogl0(x2> 

y21  =  alogi0(y2> 

x31  =  aloglO ( x3) 

y31  *  alogl0<y3) 

x41  =  alogl0(x4) 

y41  *  alogl0(y4) 

a  =  <yll-y21)/(xil-x21) 

b  =  yll  -  a#xil 

b  =  10 ,0**b 

c  =  <<y31-y41)*<x21~x41)~(x31-x41>*(y21-y41>>/ 

1  <<<x31*x31-x41*x41>*(x21-x41>>-< <x31-x41>* 

2  <x21*x21-x41*x41>>> 

d  *  <y21-y41-c*(x21*x21-x41*x41>)/<x21-x41> 
e  =  <y41-(c*x41*x41>-d*x41> 

410  if ( xlan . ge . 2 . 048) go  to  420 

c  xc  jjc  i  i  xc  x^  *  4 1  i  x^  i  x^  x^  x^  x^  x^  i  i  x^  i  i  i|t  i  i  i  xt  x^  x^  1 1 1  x  i  t  xt  x  i  i  x^  x^  x  x^  x^  3jc  x  x  x  x  x  x  x  x  x  x  x 

c*  Calculates  inpulse  if  in  parabolic  portion 

f  =  alogiO(xlan) 
g  =  <c*f*f >+<d*f >+e 
xinp  =  r ange* C 10 . 0**g ) 

c*  Calculates  inpulse  if  in  linear  portion 
goto  430 

420  xinp  =  range*b*( ( xlan)XXa> 

430  continue 

c  #  #  *  *  *  *  *  *  *  *  #  *  **  **  *  *  *  *  #*  #  *  *  *  *  *  *  #**  #  *  *  #  *  #  *  *  #  *  *  **  #  *  *  *  #  *  *  *  *  *  * 
c#  The  following  equations  are  fron  an  article  by 

c#  L.  Giglio-Tos  and  T.E.  Linnenbr ink ,  "Air  Blast 

c*  Pressure  Measurenent  Systens  and  Techniques", 
c*  MINUTES,  15th  DODESB  Seninar,  Septenber  1973, 
c*  p.  1400.  Although  the  equations  in  the  reference 
c#  have  considerable  overlap,  Lanbda  =  71.725155  was 

c*  choosen  as  the  change  over  point  because  this  is 

c#  the  point  where  the  two  equations  are  equal  (to 
c*  9  decinals) . 

c*  The  first  equation  calculates  the  peak  overpressure 
c*  <pn)  in  the  range  0 . 5<pn<71 . 725155 

c ******************************************************** 

if(xlan.gt .71 . 725155) go  to  440 

xlnpn  ■  7. 0452041-1 . 6277561*<alog< xlan> >-0 .27399088* 

1  ( ( a  1 og ( xlan) >**2) 

2  -0 . 0659731 36* ( (alog(xlan) )**3)+0 . 0065412563* 

3  < (a  1 og < x lan ) )**4) 


T 
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4  +0 . 048236359* < (alog ( xlan) )*#5)-0 . 020072553* 

5  ((alog(xlan) >**6> 

6  +0 . 0 03019 0449* ( (alog(xlan) >**7>-0 . 00015984026* 

7  ( ( al og ( xlan) >**8) 
pn  =  exp(xlnpn) 

go  to  450 

c******************************************************** 
c*  Calculates  peak  overpressure  (pn)  in  the  range 
c*  71 . 725155 (pn< 1000 . 0 

c ******************************************************** 

440  pn  *  226. 61762/( (xlan)**(i .4065913)) 
c******************** ************************************ 
c*  Calculates  the  tine  of  the  positive  phase  of  the 
c*  blast  wave  (xto)  fron  an  analytical  integral 

c*  solution  of  the  ideal  blast  wave. 

c  *********************  ik*******************]*}*******]******* 

450  xto  *  0 . 001*xinp*2 . 7 1828 1828/pm 
c ******************************************************** 
c*  Operator  inputs  length  of  ideal  blast  wave,  should 
c*  natch  length  of  actual  blast  wave 

display  'nunber  of  points  in  ideal  blast  wave  ' 
accept  len 

c*  Operator  inputs  the  sanple  interval  of  the  actual 
c*  data,  this  correlates  the  ideal  and  actual  data. 

display  'input  tine  per  point  ' 
accept  tine 
t  =0.0 

c ******************************************************** 
c*  Operator  chooses  the  desired  %  of  peak  overpressure 
c#  This  allows  hin  to  natch  the  low-pass  filter  cutoff 
c*  used  on  the  analog  data. 

c ******************************************************** 
display ' Input  desired  X  of  peak  (95-100)' 
accept  percent 
if(percent  eq,100)go  to  149 
if (percen' . eq . 99)a=410 ,87*pii/xto 
if (percent . eq . 98) a  =  178 . 58*pii/x  to 
if(percent . eq . 97)a=10B . 273*pi i/x t o 
if (percent .eq . 96) a =75 ,339*pii/xto 
if (percent . eq . 95) a*56 . 5632*pii/xt o 
ato  *  a*xto/(a*xto-l . 0 > 

c ******************************************************** 
c*  Calculates  and  files  P(t>  if  95-99%  is  chosen 
c******************************************************** 
do  500  j=i,len 
b=  -t/xto 
c=  exp<-a*t> 
d=exp (b ) 

x(j)  =  (ato*((l  0+b)*d-c)+(d-c)*(ato/(a*xto-i .0>)>*pn 
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y(j)  *  0.0 
t  =  t+tinc 
write(25>  x  <  j ) 

500  continue 

go  to  151 

^  ||^  ||^  j|^  j|^  j||  j|^  j|^  ||^  j|^  |||  ^  j|^  ^  ^  |||  ^  j||  ^  ^  ^  ^  ^  ^  ^  j|^  |||  ^  ^  ^  ^  ^  j||  ^ 

c#  Calculates  and  files  P(t>  if  100%  is  chosen 
c******************************************************** 

149  do  ISO  j  =  i,len 

x<j)  =  <p«/<exp < t/xto > > >*< i-t/xto ) 
y(j)  =0.0 
t  =  t  +  tine 
wr i te<25) x  <  j  > 

150  continue 

c  ****  *  ****  *  **********  ****  *****  *  ******  **  *******  *  ****  *****  Dc 
c#  Zeros  renainder  of  array  for  FFT 

c ******************************************************** 

151  do  155  j=  ( len  +  i>,  n 

x<j)  =  0.0 

y(j)  =  0.0 

155  write<25>  x  <  j  > 

endfile  25 

c ******************************************************** 
c*  Takes  FFT  of  ideal  blast  wave 

c ******************************************************** 

call  nuf f t ( x , y , i2p , i ) 
do  160  f  =  i,n 
wr i te(26 ) x  <  f  > 
wr ite<27) y  <f ) 

160  continue 

endfile  26 
endfile  27 

c*  Reads  in  actual  blast  wave, 
c ************************************ 
do  170  it  =  i,n 
r ead  <28) x  < i t ) 
y < it )  «  0.0 

170  continue 

c »***#*:  ***************************** 
c*  takes  FFT  of  actual  blast  wave 
c ************************************ 
call  nuff t< x , y , i2p  ,  i ) 
rewind  26 
rewind  27 

c ******************************************************** 
c*  On  a  point  by  point  basis,  coMputes  Magnitude  of 

c*  actual  FFT  and  then  the  Magnitude  of  the  ideal  FFT 

c*  and  coMpares  the  two.  If  ideal  is  larger  than  the 
c*  actual,  then  the  actual  value  is  used.  If  actual  is 
c*  larger  than  ideal  by  2.5X,  then  the  ideal  is  used  as 
c*  the  frequency  doMain  value  of  the  filtered  data. 
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do  190  f=i,n 

bnag  *  sqr  t  <  <  x ( f )##2>  +  ( y (f >*#2>  > 

read (26)pr 

read(27>pi 

pwag  *  sqr t < <pr*#2)+(pi##2> > 
if (bnag . It . <1 . 025#pnag> >go  to  180 
x  <f )  =  pr 
y  <  f )  ■  pi 
180  continue 
190  continue 

c*  Transforns  data  back  to  tine  donain. 
call  nuf f t < x , y , i2p ,-i > 

c#  Outputs  real  and  imaginary  parts  of  data  to  seperate 
c#  files,  real  in  file  30  and  imaginary  in  file  31 
c * ******* * * ****** * *** * * ** ***** * * #** ******** #*##♦* # #* #*#* # 
do  210  it=l,n 
write(30) x  < it ) 
write(3i)y<it) 

210  continue 

endfile  30 
endfile  31 

c*  Conputes  Min  and  Max  of  data  for  plotting  purposes, 
c  3k  t  i  t  i  i  i  t  i  )k  iic  3k  i  )k  )k  )k  i  )k  i  k  k  5k  ik  5k  5k  5k  5k  5k  5k  3k  5k  5k  5k  3k  5k  3k  5k  5k  5k  5k  5k  5k  3k  5k  5k  5k  5k  5k  5k  5k  5k  5k  5k 

xnax  =  0.0 

xMin  •  1000000.0 

yMax  =  0.0 

yMin  =  1000000.0 

do  220  i t=  i,n/2 

if < x ( it >. gt . xnax > xnax  =  x(it) 
if(x(it).lt  XMin>XMin  =  x(it> 
if<y<it).'  ynax)ynax  *  y<it) 
if ( y ( i t ) . i v . yMin > ynin  *  y(it) 

220  continue 

write(6,230)xnax 
230  f ornat < 1 x , 39h 

1  Maxinun  of  real  part  of  filtered  data  =,fl2.5) 

write<6,240) XMin 
240  f ornat < lx ,39h 

1  HiniMUM  of  real  part  of  filtered  data  =,fl2.S> 

wr ite<6,250  >  yMax 
250  f ornat < lx ,44h 

1  MaxiMUM  of  inaginary  part  of  filtered  data  =,fi2.5> 
write<6, 260  >  yMin 
f ornat (lx, 44h 

1  Mininun  of  inaginary  part  of  filtered  data  *,fi2.5> 
stop 
end 
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^  ^  ^  ^  ||^  j|^  j||  ^  j|^  ^  ^  j|^  ^  ^||  ^  j|^  ^  j|^  ^  ^  ^  j||  ^  ^  ^  ^  ^  j|^  ^  j|| 

C*  This  is  a  subroutine  to  perform  the  Fast 

C#  Fourier  Transform  (FFT)  of  a  complex  array 

C*  where¬ 
to#  XR  =  Real  array 

C*  XI  =  Imaginary  Array 

C*  M  =  The  order  of  the  FFT  <2AM) 

C*  IDIRN  -  The  direction  of  the  Transform 

C*  i  =  forward 

C*  -i  =  reverse 

C#  The  result  returns  in  XR  and  XI 

C  *  *  *  *  *  **  *  He  *  *  *  *  *  **  *  *#  *  *  *  *  *  *  *  *  *  #  *  *  *  *  *  *  *  *  # 

SUBROUTINE  NUFFT < XR , XI ,M , IDIRN ) 

DIMENSION  XR<i>,XI(i> 

N=2**M 
P 1=3 . 1415927 
DO  20  L  =  i  ,M 
LE=2**<M+i-L> 

LEi=LE/2 
UR  =  i  .  0 
UI=0 . 0 

WR=COS (P I/FLOAT (LEi>> 

WI=-SIN < P I/FLOAT (LEi >  > 

DO  20  J=i , LEi 
DO  10  I=J,N,LE 
IP=I+LEi 
TR=XR  < I ) +XR  < IP ) 

TI=XI < I ) +XI < IP ) 

TMR=Xk  < I )-XR ( IP  > 

TMI*XI < I )-XI < IP ) 

XR ( IP ) =TMR*UR-TMI*UI 
XI<IP)=TMR*UI+TMI*UR 
XR  < I ) =TR 
10  XI < I ) =TI 

TR=UR*WR-UI*WI 
UI=UR*WI+UI*WR 
20  UR=TR 

NV2=N/2 
NMi=N- 1 
J=1 

DO  30  1=1, NM1 

IF< I  .GE.  J)  GOTO  25 

TR=XR ( J  > 

TI=XI <  J ) 

XR  < J ) =XR ( I ) 

XI(J)*XI(I) 

XR  < I >=TR 
XI ( I >=TI 

25  K*NU2 

26  IF <K  GE.  J)  GOTO  30 
J=J-K 

K=K/2 
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GOTO  26 
30  J= J+K 

IFdDIRN  EQ.  i>  GOTO  77 

DO  99  1=1, N/2 

TTR=XR <  I  )/N 

TTI*XI<I)/N 

XR ( I  >=XR  <  N-I  +  i >/N 

XI<I)=XI< N-I+i) /N 

XR<N-I+i)*TTR 

XI(N-I+i >=TTI 

99  CONTINUE 
TTTR  =  XR  <  N ) 

TTTI*XI <  N  > 

DO  777  I=N ,2,-i 
XR ( I  )=XR  < I-i > 

XI < I )=XI ( I-i > 

777  CONTINUE 

XR(i)=TTTR 
XI<i)=TTTI 
77  RETURN 

END 

C******************************************************* 
C*  THIS  SUBROUTINE  SHIFTS  THE  ARRAY  TO  THE  RIGHT 

C*  N/2  SAMPLES  IN  A  CIRCULAR  SHIFT 

C******************************************************* 

SUBROUTINE  SHIFT<X,LS,NS> 

DIMENSION  X(LS> 

DO  200  J  =  i,NS 
Y  *  X(i) 

DO  100  K  *  2 ,LS 

100  X(K-i)  *  X ( K ) 

200  X(LS)  =  Y 

RETURN 

END 
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Program  SMOOTH 

c******************************************************* 
c*  This  program  was  developed  as  part  of  an  Electrical 
c#  Engineering  Thesis  by  Theodore  J  Moody,  entitled, 
c *  "Instrumenting  an  Explosive  Test  Arena,"  at  the 

c*  University  of  Utah,  May  1982. 

c******************************************************* 
♦control  nolist 
♦control  file  =  25-35 

integer  outfile 

dimension  psi(500,2>,  psit(2,500),  x(500,i>, 

1  psi2(500,2>,  thetax(2,i>,  thetay(2,l), 

2  tm2(2 ,500 ) , tmi (2 ,2) , y  <500 , i ) 
display  'input  file  number  for  Blastwave  ' 
accept  infile 

display  'input  file  number  for  filtered  wave  ' 
accept  outfile 

c***************************************************«**** 
c*  Operator  inputs  size  of  sliding  window 

10  display 

1  '  input  number  of  samples  to  "average"  (EVEN)  ' 

accept  nsan 
nsam2  =  nsam/2 

display  'number  of  samples  in  data  ' 
accept  ndata 

if (ndata . le . nsam)go  to  10 

call  a ve( inf ile, outfile ,nsam,nsam2,ndata) 
stop 
end 
c 
c 

C ******************************************************** 

c*  The  main  part  of  this  program  is  contained  in  the 
c*  following  subroutine  to  allow  for  the  proper 
c*  dimensioning  of  the  arrays. 

c*************************** ***************************** 
subroutine  ave< infile, outfile ,nsam,nsam2,nda ta ) 
integer  outfile 

dimension  psi<nsam,2>,  psit (2,nsam> ,  x(nsam,i>, 

1  psi2(nsam,2) ,  thetax<2,i),  thetay(2,l), 

2  tm2<2,nsam> , tmi<2,2> ,y <nsam,i ) 
c****»******«*********************«****»*****«*t********* 
c*  Builds  initial  psi  and  X  arrays 

c***«******«**«***************«*******************»*»**** 
do  100  j  *  l,nsam 
psi(j,l)  ■  float(j) 
psi< j,2>  ■  1 .0 
100  read<infile)x( j,l> 

c******************************************«******«****** 
c*  computes  least  square  of  the  first  section  of  data 
c4t******$**f*****f****ik**«$*»*t**»**«*$«**f*«*t****t**ttt 


call  xleast < >. ,psi ,a ,b ,nsa«) 
k  *  i 


c******************************************************** 
c*  Outputs  the  first  half  window  length  of  smoothed  data 

do  200  j  =  i,nsan2 
p  *  a#psi< j>i)  +  b 
display  k,p 
write( outf ile)p 
200  k  ■  k  +  i 
k  *  k  -  i 

c*  Starts  sliding  window  over  the  data  one  point  at  a 

c#  tine  conputing  a  new  least  square  for  each  new  data 

c*  point  and  outputs  a  point  representing  the  center  of 
c*  the  window. 

c******************************************************** 
do  iOOO  j=<nsan  +i)>nd'ta 
k  *  k  +  i 

call  shif t < x ,nsa«, i , i ) 
call  shift<psi>nsan>2>i) 
psi(nsan,i)  -  f loat ( k+nsan2) 
read ( inf i le> end-20 00  >  x (nsan , i  ) 
call  xleast < x >psi >a ,b >nsan> 
p  =  a*psi < nsan2> i >  +  b 
display  k,p 
write(outfile)p 
i 0  00  continue 

2000  do  2100  i  =  l,(nsan2) 

k  *  k  +  i 

^  ^  ^  A  ^  d.  ^  X  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  da  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  iL  a>  d>  ^  ^  da  ^  da 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ™  ^  ^  ^  ^  ^  ^  ^  v  ^  * 

c*  Conputes  and  outputs  the  last  half  window  length  of 

c#  snoothed  data  using  the  nost  recent  least  square 

c*  estinate. 

c #*#*#*##*#*#***##****#****#**#**#**#***#*#####*******#*# 
p  *  attpsi ( nsan2+i , i )  ♦  b 
display  k,p 

2100  write< outfile)p 

endfile  outfile 
ret  urn 
end 
c 
c 

C********«*********************************************** 

c*  the  following  subroutine  shifts  the  first  colunn  of 
c*  the  natrix  a,  no  places  up 

c******************************************************** 

subroutine  shift(a,n,n,no> 
dinension  a(n,n) 
do  100  j*i,<n-no) 
a< j,i>  ■  a< j+no,l> 
continue 
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return 

end 

c 

c 

c 

c 

c******************************************************** 
c  This  subroutine  coMputes  the  least  square 
c  using  a  standard  Least  Squares  technique, 
e******************************************************** 
subroutine  xleast (x ,psi ,a ,b >n) 

diwension  psi(n,2),  psit(2,n),  x(n,i),  tMi(2,2>, 
i  thetax(2,i) ,thetay(2,i) ,  y(n,i> , tn2(2,n) 

c ******************************************************** 
c  The  following  is  the  he^rt  of  the  calculation 
c  First  a  square  Matrix  is  forned  with 

c  psit  x  psi  =  t m1 

c  then  inverted  and  Multiplied  by  psit  the  result 

c  is  in  tM2.  Then  tM2  x  (the  vector  X)  =  thetax  (x 

c  conponent  of  the  solution) 

c 

c  thetax  ■  (psit  x  psi)A-i  x  psit  x  X 

c 

c  For  a  full  explaination  see  text  by  G.F.  Franklin 

c  and  J.D.  Powell,  DIGITAL  CONTROL  OF  DYNAMIC  SYSTEMS, 

c  (Addison-Wesley  Publishing  Co.,  Massachusetts, i98i) , 

c  pp  .  216-220. 

c ******************************************************** 
call  xt (psi ,psit ,n ,2) 
call  xM(psit ,psi , tMi ,2,n ,2, i . 0 > 
call  xinvert(2,tMi,ierror> 
call  XM(tMi,psit,tM2,2,2,n,i.0> 
call  xm( tn2, x , thetax ,2,n , 1 , i . 0 ) 
a  *  thetax(i,i) 
b  *  thetax ( 2,1) 
r eturn 
end 

c******************************************************** 

c 

c  TRANSPOSES  A  MATRIX 

c******************************************************** 

Subroutine  xt(a,b,M,n) 
dinension  a(M,n ) ,b (n ,m) 
do  100  i»l,n 

do  100  j  *  1,m 
b(i, j)  *  a( j,i> 

100  continue 
return 
end 
c 
c 
c 
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c******************************************************** 

c 

c  MATRIX  MULTIPLY 
C 

C******************************************************** 

c 

Subroutine  xm (a, b,c,na,na,nb, scaler) 
dinension  a(na,na),b(na,nb’/,c(na,nb) 
do  200  i  «=  i , nb 
do  200  j  *  i,na 
c<  j  ,  i )  *  0.0 
do  100  k  *  i , na 

c(j,i>  ■  c(j,i)  +a < j , k )»b < k , i ) 

100  continue 

c<j,i)  *  sea ler#c < j , i ) 

200  continue 
return 
end 
c 

c******************************************************** 
c*  The  following  exits  the  progran  if  the  Matrix  being 
c#  inverted  is  singular. 

Subroutine  exitinverse 
wr  it:><6, 100 ) 

100  f orna t < ix ,3ihMatr i x  is  singular,  check  input) 

urite<6,200) 

200  fornatdx,  12hEXIT  PROGRAM) 
stop 
end 
c 

c ******************************************************** 
c*  COMPUTES  THE  INVERSE  OF  A  MATRIX 

Subroutine  xinver t (n ,a , ierror > 
dinension  a(n,n> 
if<n.gt.2)go  to  1000 
tenpaa  <2 ,2) 
a(2,2)  *  ad ,  i  > 
a (1,1)  =  temp 
a ( i , 2)  *  -a< i ,2) 
a ( 2, i )  *  -a ( 2 , i ) 

tenp  *  ad,i)*a<2,2)  -  ad,2)*a(2,l) 
if <abs<  tenp) . It . lE-7)call  exitinverse 
do  100  i  -1,2 

do  100  j  =  i  ,2 

100  a(i,j)  «  a(i,j)/tenp 

go  to  3000 

c******************************************************** 
c*  HP3000/33  intrinsic  inversion  routine 
c  ****************************************************<•*** 

1000  call  invert <n,a d , i) , ierror ) 


87 


display  ierror 

if < ierror .lt.0.5)call  sx it in verse 
3000  return 
end 


APPENDIX  F 


PRINTOUT  OF  RANGE-BANDWIDTH 
PRODUCTS 


I 

\ 

i 

. . — .  ..  ■■■  ■  -  ■  ■■—  —n,  i  „  .11..  - - - 


Laabda 

Range-Bandwidth 
Product  (ft-Hi) 

Peak 

Laabda 

Range-Bandwidth 
Product  (ft-He) 

Peak 

. 519397 

144074. 

.900000 

3.47136 

90234.8 

.900000 

S39Sii 

152819. 

.900000 

3 . 60S79 

88711 .9 

.899999 

. 560402 

161465. 

.900000 

3.74542 

87199.6 

.900000 

.582103 

169948. 

.900000 

3.89045 

85703 . 0 

.899999 

.604644 

178200 . 

.899999 

4.04111 

84226  8 

.900000 

.628058 

186154. 

.900000 

4 . 19759 

82775.7 

.899999 

. 652379 

193742 . 

.899999 

4.36014 

81353.8 

.900000 

.677641 

200894. 

.900000 

4.52898 

79965 . 0 

.899999 

.703882 

207544. 

.900000 

4.70436 

78612.7 

.900000 

.731139 

213626. 

.900000 

4 . 88653 

77300.5 

.900000 

. 759451 

219078 . 

.900000 

5.07575 

76030.9 

.900000 

. 788860 

223841 . 

.899999 

5.27230 

74806.9 

.899999 

.819408 

227863 . 

.  900000 

5.47646 

73630 . 8 

.899999 

.851138 

231095. 

.900000 

5 . 68853 

72504.6 

.900000 

. 884097 

233501. 

.900000 

5. 90881 

71430.3 

. 899999 

. 918333 

23S048. 

.900000 

6 . 13762 

70409.5 

.899999 

. 9S3B94 

235713. 

.900000 

6 . 37S29 

69443 . 5 

.899999 

.990832 

235486 . 

.900000 

6.62217 

68533.8 

.899999 

1.02920 

234362 . 

.899999 

6.87860 

67681 . 2 

.900000 

1.06905 

232350 

.900000 

7 . 14497 

66886 . 8 

.900000 

1 .11045 

229469 . 

.899999 

7.42164 

66151.2 

.899999 

1 . 15345 

22S745 

.899999 

7.70904 

65475.3 

.900000 

1 .19812 

221218 

.900000 

8  00756 

64859 . 4 

.899999 

1.24451 

215935. 

. 899999 

8  31764 

64304.1 

.900000 

1.29271 

209951 . 

.900000 

8 . 63973 

63809  7 

. 899999 

1.34277 

203331 . 

899999 

8 . 97429 

63376 . 6 

.899999 

1 . 39476 

196142. 

. 899999 

9.32181 

63005.0 

.900000 

1  44877 

188461 

.900000 

9 . 68279 

62695.2 

.900000 

1  S0487 

180363. 

.900000 

10.0577 

62447  5 

.900000 

1 . 56315 

171931 . 

. 899999 

10.4472 

62261  8 

.899999 

1 . 62368 

163246  . 

.900000 

10.8518 

62138  6 

. 899999 

1.68655 

1S43B7 . 

. 899999 

11.2720 

62077.9 

.900000 

1.75186 

145433. 

.899999 

11.7085 

62079  9 

.900000 

1  81970 

136460. 

.900000 

12.1619 

62145.0 

.900000 

1 .89017 

127S40 . 

. 899999 

12.6328 

62273  0 

.900000 

1 .96336 

118738 

.900000 

13.1220 

62464.3 

.899999 

2 . 03939 

110115. 

. 899999 

13  6301 

62719.3 

.900000 

2.11836 

108255 

.899999 

14.1579 

63038.1 

.900000 

2.20039 

107153. 

.900000 

14.7062 

63420.7 

.899999 

2.28560 

105979. 

. 899999 

15.2757 

63867 . 6 

.900000 

2.37411 

104741 . 

.899999 

15.8672 

64379 . 0 

.899999 

2  46604 

103445. 

.899999 

16.4816 

64955 . 1 

. 899999 

2.S61S3 

102097. 

. 899999 

17.1199 

65596 . 1 

.900000 

2.66073 

100704. 

.900000 

17.7828 

66302  3 

. 899999 

2 . 76376 

99272 . 5 

. 899999 

18.4714 

67074 . 0 

.900000 

2.87078 

97808.7 

.899999 

19.1867 

67911  2 

.900000 

2 . 98195 

96319.4 

. 899999 

19.9297 

68814  0 

.900000 

3 . 09742 

94811.1 

.900000 

20.7014 

69782 . 8 

.900000 

3.21736 

93290 . 2 

.900000 

21.5031 

70817.0 

.900000 

3.34195 

91762  6 

.900000 

22 . 3357 

71917.1 

.900000 
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Laabda 

Range-Bandwidth 
Product  (ft-Hs) 

Peak 

Laabda 

Range-Bandwidth 
Product  (ft-Hz) 

23.2007 

73083.0 

.899999 

155  060 

183820 

24.0991 

74314.4 

.900000 

161.065 

186794. 

25 . 0323 

75610.7 

.899999 

167.302 

189616. 

26.0016 

76972 . 3 

.900000 

173.780 

192887 . 

27.0085 

78397 . 8 

.899999 

180.510 

196007. 

28.0544 

79886 . 6 

. 899999 

187.500 

199178. 

29.1407 

81438.8 

.899999 

194.760 

202400 . 

30.2692 

83052.5 

.900000 

202.302 

205675 . 

31.4413 

84728 . 2 

.899999 

210.136 

209002 

32 . 6588 

86462.5 

.899999 

218.273 

212383 

33  9235 

88255 . 5 

.900000 

226.725 

215819. 

35.2371 

90106.6 

.899999 

235.505 

219311  . 

36.6016 

92012.8 

.899999 

244.625 

222859 . 

38.0190 

93971.5 

.899999 

254.097 

226464 . 

39.4912 

95982 . 4 

.899999 

263 . 937 

230128. 

41.0205 

98044.3 

.900000 

274.158 

233851 . 

42.6089 

100153. 

.900000 

284  774 

237634 

44.2589 

102306. 

.900000 

295.802 

241479. 

45 . 9728 

104505. 

. 899999 

307.256 

245385 . 

47.7530 

106744. 

.900000 

319.154 

249355 . 

49  6022 

109022. 

.899999 

331.513 

253389 

51.5229 

111336. 

.900000 

344.350 

257489 

S3. 5181 

113685. 

.899999 

357.685 

261654. 

55.5905 

116064. 

. 899999 

371.535 

265887 . 

57.7431 

118480 . 

.900000 

385 . 923 

270188 

59  9791 

120915 

.900000 

400.867 

274559 . 

62  3018 

123382. 

. 899999 

416.390 

279001 

64.7144 

125869. 

.900000 

432.514 

283515 

67  2203 

128387 . 

. 899999 

449.262 

288102. 

69.8233 

130920. 

. 899999 

466  659 

292762 . 

72  5271 

133351 . 

099999 

484.731 

297499 . 

75  3357 

135509. 

.899999 

503.501 

302311 . 

78.2529 

137701 

.900000 

522.998 

307202. 

81  .2831 

139929. 

.900000 

543.250 

312172. 

64.4307 

142193. 

.900000 

564 . 287 

317222. 

87.7002 

144493. 

.900000 

586.138 

322354 . 

91.0963 

146831 . 

.900000 

608.835 

327569 

94 . 6238 

149206. 

. 899999 

632.412 

332869 . 

98.2880 

151620 . 

.899999 

656.902 

338254 . 

102.094 

154073. 

.900000 

682.339 

343726. 

106.048 

156565. 

.900000 

708.761 

349287. 

110.154 

159098. 

.900000 

736.208 

354938 . 

114.420 

161672. 

.900000 

764.717 

360681 . 

118  850 

164287 . 

.900000 

794.330 

366515. 

123.453 

166945. 

. 899999 

825.089 

37244S . 

128.233 

169646. 

.900000 

857.039 

378470 . 

133.199 

172391 . 

. 899999 

890.227 

384593 

138  357 

175180 . 

.900000 

924 . 699 

390816. 

143.714 

178014. 

.900000 

960  508 

397137 

149.280 

180894. 

.900000 

997.702 

403562. 

Peak 

099999 
.899999 
.900000 
.900000 
.900000 
.  899999 
.900000 
.900000 
.899999 
.900000 
.899999 
.899999 
.900000 
.900000 
.899999 
.899999 
.900000 
.  899999 
.899999 
.900000 
.900000 
.899999 
.899999 
.899999 
.900000 
.899999 
.900000 
.900000 
.900000 
.899999 
.900000 
. 899999 
.900000 
.900000 
.899999 
.900000 
.900000 
.899999 
.899999 
.900000 
.899999 
.900000 
.900000 
.900000 
.900000 
.900000 
.899999 
.899999 
.899999 
.900000 
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Lambda 

Range-Bandwidth 
Product  <ft-Hz) 

Peak 

Lambda 

Range-Bandwidth 
Product  (ft-Ha) 

Peak 

. S19397 

36S462 . 

.950000 

3.47136 

228893 

.950000 

. S39S11 

387647 . 

.950000 

3 . 60579 

225030 . 

.950000 

.560402 

409579 . 

.950000 

3 . 74542 

221194. 

.950000 

.582103 

431095 

.950000 

3.69045 

217397 

.950000 

.  604644 

452029 . 

.950000 

4.04111 

213653 

.950000 

. 628056 

472205 . 

.950000 

4.19759 

209972. 

.950000 

. 652379 

491452. 

.950000 

4.36014 

206365 . 

.950000 

. 677641 

509596 

.950000 

4 . 52898 

202842. 

.950000 

.703862 

526465 . 

.950000 

4.70436 

199412. 

.950000 

.731139 

541692. 

.950000 

4.88653 

196083. 

.950000 

. 7574S1 

555721 . 

.950000 

5 . 07575 

192863. 

.950000 

.788860 

567803 . 

.950000 

5 . 27230 

189758. 

.950000 

.819408 

578005. 

.9S0000 

5 . 47646 

186775 . 

.950000 

.851138 

566205. 

.950000 

5 . 68853 

183912 

.950000 

.884097 

592307. 

.950000 

5 . 90881 

181193. 

.950000 

.918333 

596231 . 

.950000 

6.13762 

178603. 

950000 

. 9S3894 

597919  . 

.950000 

6 . 37529 

176153 

.950000 

. 990832 

597343 

.950000 

6.62217 

173845 

.950000 

1.02920 

594492 

. 9S0  000 

6 . 87860 

171683 . 

950000 

1 . 0690S 

589388 

.950000 

7.14497 

169668. 

950000 

1  11045 

582079 . 

.950000 

7.42164 

167802. 

.950000 

1 . 1S345 

572634 

.950000 

7.70904 

166087. 

.950000 

1.19812 

561149 

.950000 

8.00756 

164525 

,950000 

1  244S1 

547748 . 

.950000 

8.31764 

163116. 

.950000 

1.29271 

532570 . 

.950000 

8 . 63973 

161862. 

.950000 

1 .34277 

515776 

.950000 

8.97429 

160764 

.950000 

1.39476 

497S41 . 

. 950000 

9.321B1 

159821 

.950000 

1 . 44877 

478056 . 

.950000 

9.68279 

159035 . 

.950000 

1 . 50487 

457516. 

.950000 

1 0 . 0577 

158407. 

.950000 

i  56315 

436128. 

.950000 

10.4472 

157936 

.950000 

1  62368 

414096 

.950000 

10.8518 

157623 

.950000 

1 . 68655 

391623. 

.950000 

11.2720 

157469 . 

.950000 

1  7S186 

368912. 

.950000 

11 .7085 

157474 

.950000 

1.81970 

346150. 

.950000 

12  1619 

157639 

.950000 

i  89017 

323522 

.950000 

12  6328 

157964 

.950000 

1 . 96336 

301195. 

.950000 

13  1220 

158449 . 

.950000 

2 . 03939 

279321 . 

.950000 

13  6301 

159096 

.950000 

2.11836 

274604 . 

.950000 

14.1579 

159905 

.950000 

2.20039 

271807. 

.950000 

14.7062 

160875 

.950000 

2 . 28560 

268831 . 

.950000 

15.2757 

162009 

.950000 

2.37411 

265691 . 

.950000 

15.8672 

163306. 

950000 

2 . 46604 

262403. 

.950000 

16.4816 

164767 . 

.950000 

2.56153 

258984 . 

.950000 

17.1199 

166393 

.9S0000 

2 . 66073 

255450 . 

.950000 

17.7828 

168185. 

.950000 

2.76376 

251818. 

.950000 

18.4714 

170143. 

.950000 

2 . 87078 

248105. 

.950000 

19.1867 

172266 

.950000 

2.98195 

244327 . 

.950000 

19.9297 

174556 . 

.950000 

3 . 09742 

240501. 

.950000 

20.7014 

177014. 

.950000 

3.21736 

236643. 

.950000 

21  .5031 

179637. 

.950000 

3  3419S 

232768. 

.950000 

22.3357 

182428 . 

.950000 

Lambda 

Range- Bandwidth 
Product  (ft-Bz) 

Peak 

Lambda 

Range-Bandwidth 
Product  (£t-Hz) 

Peak 

23.2007 

185385. 

.950000 

155.060 

466284 . 

.950000 

24.0991 

188509 . 

.950000 

161.065 

473828. 

.950000 

25.0323 

191797. 

.950000 

167.302 

481494 

.950000 

26.0016 

195251 . 

.950000 

173.780 

489283 . 

.950000 

27  0085 

198867. 

.950000 

180.510 

497199. 

.950000 

28.0544 

202643 

.950000 

187.500 

505243. 

950000 

29.1407 

206581 

.950000 

194.760 

513416 

.950000 

30.2692 

210674 

.950000 

202.302 

521722. 

.950000 

31.4413 

214925. 

. 950000 

210.136 

530163 

.950000 

32  6588 

219324 

.950000 

218  273 

538739 

.950000 

33  9235 

223872. 

.950000 

226.725 

547455 

950000 

35.2371 

228568 . 

.950000 

235  505 

556313 

.950000 

36  6016 

233403 

.950000 

244.625 

565311 

.950000 

38.0190 

238372 

.950000 

254.097 

574457 

950000 

39.4912 

243473 

.950000 

263.937 

583751 . 

950000 

41 . 0205 

248703. 

. 950000 

274. 158 

593195. 

.950000 

42  6089 

254051 . 

.950000 

284  774 

602792 

9SOOOO 

44  2589 

259513. 

. 950000 

295.802 

612544 

.950000 

4S  9728 

26S092. 

.950000 

307.256 

622452 

.950000 

47 .7530 

270770 

.950000 

319  154 

632523 . 

950000 

49  6022 

276548 

.950000 

331 .513 

642756 

.950000 

SI  .5229 

282420 

.950000 

344  350 

6S315S 

9SOOOO 

53.5181 

288377 

.950000 

357.685 

663721 

950000 

SS  S90S 

294413 

.950000 

371 .535 

674459. 

.950000 

57.7431 

300541 

.950000 

385  923 

685370 

.950000 

59.9791 

306719 

.950000 

400  867 

696457 . 

950000 

62.3018 

312976 

.950000 

416.390 

707725 

.950000 

64 .7144 

319283 

.950000 

432.514 

719175. 

950000 

67  2203 

325672 

.950000 

449  262 

730810 

.950000 

69  8233 

332097 

.950000 

466  659 

742632. 

.950000 

72  5271 

338264 

.950000 

484.731 

754646 

.950000 

75.3357 

343737 

.950000 

503  501 

766853 

.950000 

78  2529 

349298 

.950000 

522.998 

77V2(.x 

.950000 

81  .2831 

354949 . 

.950000 

543.250 

791867. 

.950000 

84  4307 

360691  . 

.950000 

564.287 

804677. 

.950000 

87.7002 

366527 . 

.950000 

586  138 

817696 

.950000 

91  0963 

372456 

.950000 

608  835 

B30923 . 

.950000 

94  6230 

378481 . 

.950000 

632  412 

844367. 

.950000 

98  2880 

384605 

.950000 

656  902 

858029 

950000 

102. 094 

390826 

.950000 

682  339 

871909 

.950000 

106.048 

397150  . 

950000 

708.761 

886015. 

.950000 

110.154 

403575. 

.950000 

736.208 

900349. 

.950000 

114.420 

410104 . 

.950000 

764.717 

914916. 

.950000 

118  850 

416737 

.950000 

794.330 

929716 

.950000 

123  453 

423480 

.950000 

825.089 

944759 

.950000 

128  233 

430331  . 

.  950000 

BS7. 039 

960041 

.950000 

133  199 

437293 

.950000 

890.227 

975573 . 

.950000 

138  357 

444367 . 

.950000 

924.699 

991358 . 

.950000 

143.714 

451556 

.950000 

960.508 

.  100739E+07 

.950000 

149.280 

458862 . 

.950000 

997.702 

.  102369E+07 

950000 

Laabda 

Range-Bandwidth 
Product  (ft- Bs) 

Peak 

Laabda 

Range- Bandwidth 
Product  (ft-Mz) 

Paak 

.519397 

486775 . 

.960000 

3.47136 

304872. 

.960000 

.539511 

516324. 

.960001 

3.60S79 

299727. 

.960000 

.550402 

545536 

.960001 

3 . 74542 

294618 . 

.960000 

.582103 

574195. 

.960000 

3.89045 

269561 . 

.960000 

.504544 

602077. 

.960000 

4.04111 

284574 . 

.960000 

. 528058 

628951 . 

.960000 

4 . 19759 

279671 . 

.960001 

.552379 

654587 . 

.960000 

4.36014 

274867 . 

.960001 

.677541 

678753 

.960000 

4 . 52898 

270174. 

.960000 

.703882 

701221 . 

.960000 

4.70436 

265606. 

.960001 

.731139 

721770. 

. 960000 

4 . 88653 

261172. 

.960000 

.759451 

740189 

.960000 

5.07575 

256883 . 

.960000 

. 788850 

756281 . 

. 960000 

5.27230 

252747 . 

.960001 

.819408 

769870. 

.960000 

5.47646 

248773 . 

.960000 

.851138 

780792 

.960000 

5 . 68853 

244968. 

.960001 

.884097 

788920 . 

.960000 

5.90881 

241339. 

.960000 

.918333 

794147. 

.960000 

6.13762 

237890 

.960001 

9S3894 

796395 . 

.960000 

6 . 37S29 

234626 . 

.960001 

.990832 

795627 . 

.960000 

6.62217 

231552 . 

.960000 

1.02920 

791830 . 

.960000 

6.67860 

228672 . 

.960001 

1 .06905 

785032 

.960000 

7 . 14497 

225988 

.960000 

1.11045 

775296 . 

.960000 

7.42164 

223502. 

.960000 

1 . 15345 

762716 

.960000 

7.70904 

221219. 

.960000 

1.19812 

747418. 

.960001 

8.00756 

219138 

.960000 

1 .24451 

729570 . 

.960000 

B. 31764 

217262. 

.960001 

1  29271 

709353. 

.960001 

B . 63973 

215591 . 

.960000 

1 . 34277 

686985 . 

.960000 

B  97429 

214128. 

.960000 

1  39476 

662697 . 

.960001 

9  32181 

212872 

.960000 

1 . 44877 

636744 

.960001 

9 . 68279 

211826. 

.960000 

1.50407 

609385. 

.960000 

1 0 . 0577 

210989 

.960000 

1.56315 

580898 

.960000 

10.4472 

210361. 

.960001 

1 . 62368 

551552. 

.960000 

10.8518 

209945. 

.960001 

1 . 68655 

521620. 

.960000 

11.2720 

209740 

.960000 

1  75186 

491369 

.960000 

11.7085 

209747. 

960001 

1.81970 

461053. 

.960000 

12.1619 

209967. 

.960000 

1  89017 

430913. 

.960000 

12.6328 

210399 

.960000 

1 . 96336 

401175. 

.960000 

13.1220 

211046. 

.960000 

2.03939 

372040 . 

.960001 

13.6301 

211907. 

.960000 

2.11836 

365757 

.960000 

14.1579 

212984  . 

.960000 

2.20Q39 

362032 . 

.960000 

14.7062 

214277. 

.960000 

2  2B560 

358068 

.960000 

15  2757 

215787 . 

.960000 

2.37411 

353885 . 

.960000 

15.8672 

217515 

.960000 

2.46604 

349507. 

.960000 

16 . 4816 

219461 

.960000 

2  56153 

344952 . 

.960000 

17.1199 

221627 

.960000 

2  66073 

340245 

.960000 

17 .7828 

224013. 

.960000 

2  76376 

335408. 

.960000 

18.4714 

226620 . 

.960000 

2.87078 

330462. 

.960001 

19 . 1867 

229449. 

.960000 

2  98195 

325430 . 

.960000 

19.9297 

232499 . 

.960000 

3 . 09742 

320334 . 

.960000 

20.7014 

235772. 

.960000 

3  21736 

315196. 

.960000 

21.5031 

239266 . 

.960000 

3.34195 

310035. 

.960000 

22.3357 

242983. 

.960000 

94 


Laabda 

Range-Bandwidth 
Product  (ft-Bx) 

Peak 

Laabda 

Range-Bandwidth 
Product  < ft-Bx) 

Peak 

23.2007 

246922 . 

.960000 

155.060 

621065 

.960000 

24 .0991 

2S1083. 

.960000 

161.065 

631112 

.960000 

25.0323 

255463 . 

.960000 

167.302 

641323 

.960000 

26.0016 

260063. 

.960001 

173.780 

651698 . 

.960000 

27.0085 

264879 . 

.960001 

180.510 

662241 . 

.960001 

28.0S44 

269909 . 

.960000 

187.500 

672955 . 

.960000 

29.1407 

275154. 

.960000 

194.760 

683841 . 

.960001 

30.2692 

280606. 

.960000 

202.302 

694904 

.960000 

31.4413 

286268 . 

.960000 

210.136 

706147. 

.960000 

32 . 6588 

292127. 

.960000 

218.273 

717570. 

.960000 

33.9235 

298185. 

.960000 

226 . 725 

729179 

.960000 

35.2371 

304439. 

.960001 

235.505 

740977 . 

.960000 

36.6016 

310880. 

.960000 

244 . 625 

752963 . 

.960000 

38.0190 

317498 . 

.960000 

254.097 

765144 

.960001 

39.4912 

324292. 

.960000 

263.937 

777523 . 

.960000 

41 .0205 

331258. 

.960000 

274 . 158 

790102. 

.960001 

42.6089 

338382. 

.960000 

284.774 

802885 

.960000 

44.2589 

345657. 

.960000 

29S.802 

815874 . 

.960001 

4S . 9728 

353087. 

.960000 

307.256 

829072 . 

.960001 

47.7530 

360651 . 

.960001 

319.154 

842485 . 

.960001 

49.6022 

368347. 

960000 

331 .513 

856114 . 

.960000 

SI  .5229 

376168. 

.960001 

344.350 

869966 . 

.960000 

S3  Sl8i 

384102. 

.960000 

357 . 685 

884039. 

.960001 

SS.S90S 

392142 . 

.960000 

371.535 

898342 . 

.960000 

57.7431 

400303. 

.960000 

385 . 923 

912874. 

.960001 

59 . 9791 

408532. 

.960000 

400. B67 

927642 

.960001 

62.3018 

416866 

.960001 

416.390 

942650 . 

.960000 

64  7144 

425267. 

.960000 

432  514 

957901 

.960000 

67  2203 

433777 

.960000 

449  262 

97339B 

.960001 

69.8233 

442335 . 

.960000 

466.659 

989144 

.960001 

72  S271 

450549. 

.960000 

484.731 

. 100515E+07 

.960000 

75 . 33S7 

457839 

.960001 

503.501 

.  102141E+07 

.960000 

78  2529 

465245 

.960001 

522.998 

.  103793E+07 

.960000 

81.2831 

472772 . 

.960000 

543.250 

105472E+07 

.960001 

84.4307 

480421 . 

.960000 

564 . 287 

.  107179E+07 

.960000 

87.7002 

488193. 

.960000 

586.138 

. 108912E+07 

.960000 

91.0963 

496090 . 

.960000 

608.835 

.  110674E+07 

.960000 

94.6238 

504115. 

.960000 

632.412 

112465E+07 

.960000 

98.2880 

512272. 

.960001 

656.902 

. 11428SE+07 

.960000 

102.094 

520558. 

.960001 

682.339 

.  116133E+07 

.960001 

106.048 

528981 . 

.960000 

708.761 

118012E+07 

.960000 

110.154 

537539 . 

.960000 

736  208 

. 119921E+07 

.960000 

114.420 

546235 . 

.960000 

764.717 

. 121862E+07 

.960001 

118.850 

555071 . 

.960000 

794.330 

.  123833E+07 

.960000 

123.453 

564052. 

.960000 

825 . 089 

.  125837E+07 

.960000 

128.233 

S73177 . 

.960000 

857 . 039 

.  127872Et-07 

.960001 

133.199 

582449. 

.960000 

890.227 

. 129941E+07 

.960000 

138.357 

591872 . 

.960001 

924.699 

.  132043E+07 

.960000 

143.714 

601447. 

.960000 

960.508 

134179E*07 

.960000 

149.280 

611179. 

.960000 

997.702 

1363S0E+O7 

.960000 

95 


Lambda 

Range-Bandwidth 
Product  (ft-Bz) 

Peak 

Lambda 

Range-Bandwidth 
Product  (ft-Hz) 

Peak 

.519397 

699566 . 

.970000 

3. 

47136 

438146. 

.970000 

.539511 

742032. 

.970000 

3. 

60579 

430751. 

.970000 

.560402 

784014. 

.970000 

3. 

74542 

423408. 

.970000 

.562103 

825201 . 

.970000 

3. 

89045 

416141 . 

.970000 

.604644 

865271 . 

.970000 

4. 

04111 

408973. 

.970000 

. 628056 

903893. 

.970000 

4. 

19759 

401927. 

.970000 

.652379 

940736. 

.970000 

4. 

36014 

395023 

.970000 

. 677641 

975466 . 

.970000 

4. 

52898 

388279 . 

.970000 

.703682 

. 100776E+07 

.970000 

4. 

70436 

381713. 

.970000 

.731139 

. 103729E+07 

.970000 

4. 

88653 

37S342 . 

.970000 

.759451 

. 1 06376E+07 

.970000 

5. 

07575 

369177. 

.970000 

.788860 

. 10868BE+07 

.970000 

5. 

27230 

363234 . 

.970000 

.819408 

. 110641E+07 

.970000 

5. 

47646 

357523 . 

.970000 

.851138 

. 112211E+07 

.970000 

5. 

68853 

352055. 

.970000 

.884097 

. 113379E+07 

.970000 

5. 

90881 

346838 . 

.970000 

.918333 

. 114130E+07 

.970000 

6. 

13762 

341882. 

.970000 

. 953894 

. 1 14453E+07 

. 970000 

6. 

37529 

337191 . 

.970000 

.990832 

.  114343E+07 

.970000 

6. 

62217 

332774 . 

.970000 

1.02920 

. 113797E+07 

.970000 

6. 

87860 

328634 

.970000 

1 . 06905 

.  112820E+07 

.970000 

7. 

14497 

324777 

.970000 

1.11045 

. 111421E+07 

.970000 

7. 

42164 

321205. 

.970000 

1 . 1S345 

. 109613E+07 

.970000 

7  . 

70904 

317923. 

.970000 

1.19812 

. 107415E+07 

.970000 

8 

00756 

314932. 

.970000 

1.24451 

. 104850E+07 

.970000 

B. 

31764 

312236. 

.970000 

1.29271 

. 101944E+07 

.970000 

8. 

63973 

309836. 

.970000 

1 .34277 

987296 . 

.970000 

8. 

97429 

307733. 

.970000 

1.39476 

952391 . 

.970000 

9. 

321B1 

305928. 

.970000 

1.44877 

915093. 

.970000 

9. 

68279 

304424. 

.970000 

1.50487 

875775. 

.970000 

10 

0577 

303221 . 

.  97  C  0  0  0 

1.56315 

834833. 

.970000 

10 

.4472 

302320. 

.970000 

1.62368 

792660 . 

.970000 

10 

.8518 

301721 

.970000 

1 . 68655 

749643. 

.970000 

11 

.2720 

301427. 

.970000 

1.7S186 

706169. 

,970000 

11 

.7085 

301436 

.970000 

1.81970 

662599 

.970000 

12 

.1619 

301752. 

.970000 

1.89017 

619284. 

.970000 

12 

.6328 

302374 

.970000 

1  .96336 

576546. 

.970000 

13 

.1220 

303303. 

.970000 

2.03939 

534675 . 

.970000 

13 

.6301 

304541 

.970000 

2.11836 

525645. 

.970000 

14 

.1579 

306089. 

.970000 

2.20039 

520292. 

.970000 

14 

.7062 

307947. 

.970000 

2 . 28560 

514595. 

.970000 

15 

2757 

310117. 

.970000 

2  37411 

508584. 

.970000 

15 

.8672 

312600 . 

.970000 

2.46604 

502291 . 

.970000 

16 

.4816 

315397. 

.970000 

2.56153 

495746. 

.970000 

17 

.1199 

318509 . 

.970000 

2 . 66073 

488981 . 

.970000 

17 

.7828 

321939 . 

.970000 

2 . 76376 

482029. 

.970000 

18 

.4714 

325686 . 

.970000 

2.87078 

474922. 

.970000 

19 

.1867 

329751 . 

.970000 

2.98195 

467690 . 

.970000 

19 

.9297 

334135 

.970000 

3 . 09742 

460366. 

.970000 

20 

.7014 

338838 

.970000 

3.21736 

452982 . 

.970000 

21 

.5031 

343860 . 

.970000 

3.34195 

445564 

.970000 

22 

.3357 

349202. 

.970000 

96 


Laabda 

Range-Bandwidth 

Product 

Pul 

Laabda 

Range- Bandwidth 
Product  (ft-Hz) 

Peak 

23.2007 

354863 . 

.970000 

155.060 

B92SS9 . 

.970000 

24 . 0771 

360842 . 

.970000 

161.065 

906979 . 

.970000 

25.0323 

367137. 

.970000 

167 .302 

721674 

.970000 

26.0016 

373748 . 

.970000 

173.780 

936584 . 

.970000 

27.0085 

380670  . 

.970000 

180  510 

951736 . 

.770000 

28.0544 

387879 . 

.970000 

187.500 

967133. 

.970000 

27.1407 

39S436 

.970000 

194.760 

982778 . 

.770000 

30.2672 

403271 . 

.970000 

202.302 

978677 . 

.770000 

31.4413 

411408 . 

.970000 

210.136 

. 1014832*07 

.970000 

32.6588 

419829 . 

.970000 

218.273 

10312SE+07 

.970000 

33  7235 

428535 . 

.770000 

226 . 725 

. 104774E+07 

.970000 

35.2371 

437523 

.770000 

235.505 

. 106489E+07 

.970000 

36.6016 

446779 . 

.970000 

244.625 

108212E*07 

.970000 

38.0170 

456290 . 

.970000 

254  097 

.  109762E*07 

.970000 

37.4712 

466054. 

.970000 

263.937 

. 111741E+07 

.770000 

41.0205 

476066. 

.970000 

274  158 

.  113547E+07 

.970000 

42.6087 

486303 . 

.970000 

284.774 

.  11S386E*07 

.770000 

44.2587 

476759 . 

.970000 

295.802 

. 117253E+07 

.970000 

45 . 7728 

507437. 

.770000 

307.256 

. 117150E*0? 

.970000 

47.7530 

518307. 

.970000 

317.154 

. 121077E+07 

.970000 

47  6022 

527367 

.770000 

331  513 

. 123036E+07 

.970000 

51  S227 

540607. 

.770000 

344.350 

. 12S027E+07 

.970000 

S3  5181 

552010 

.770000 

357  685 

. 127049E*07 

.970000 

55.5705 

563564 . 

.970000 

371.535 

12710SE+07 

.970000 

57  7431 

575293 

.970000 

385  723 

131173E+07 

.970000 

59.7771 

587117. 

.970000 

400.867 

.  133316E*07 

.770000 

62.3018 

597096. 

.970000 

416.390 

. 135472E*07 

.770000 

64.7144 

611170  . 

.970000 

432  514 

. 137664E*07 

.770000 

67.2203 

623400 . 

.970000 

447.262 

. 137871E*07 

.770000 

69.8233 

635679 . 

.770000 

466.657 

142154E*07 

.770000 

72  5271 

647504 . 

.770000 

484.731 

. 144454E+07 

.770000 

75  3357 

657780 . 

.770000 

503.501 

146791E*07 

.970000 

78  2527 

668624 

.770000 

522  978 

147166E*07 

.770000 

81 .2831 

677442 

.970000 

543.250 

151577E*07 

.770000 

84.4307 

670434. 

.770000 

564 . 287 

.  154031E+07 

.770000 

87.7002 

701603. 

.970000 

586 . 138 

156523E+07 

.770000 

91 . 0963 

712953. 

.770000 

608  835 

. 157055E*07 

.770000 

94.6238 

724486. 

.970000 

632.412 

. 161628E*07 

.770000 

982880 

736208. 

.970000 

656.902 

, 164243E+07 

.770000 

102.094 

748117. 

.970000 

682 . 339 

. 166700E*07 

.770000 

106.048 

760222 . 

.770000 

708.761 

. 167600E*07 

.970000 

110.154 

772521 . 

.770000 

736.208 

. 172344E+07 

.770000 

114.420 

785017. 

.770000 

764.717 

. 175133E*07 

.770000 

118  850 

777717. 

.770000 

794.330 

.  177766E*07 

.770000 

123.453 

810623. 

.970000 

825  087 

. 180845E*07 

.970000 

128.233 

823738 . 

.770000 

857 . 037 

.  183771E+07 

.770000 

133  177 

837063 

.970000 

870.227 

. 186744E+07 

.770000 

138.357 

850605. 

.770000 

724,677 

. 187765E+07 

.770000 

143.714 

864366 . 

.770000 

960.508 

. 172835E*07 

.770000 

147.280 

878351 . 

.770000 

777.702 

. 175754E+07 

.970000 

Lambda 

Range-Bandwidth 
Product  (ft-Ss) 

Peak 

Lambda 

Ranae-Bandvi dt  h 
Product  (ft-Hs) 

Peak 

519397 

.115383E+07 

.980000 

3.47136 

722655 . 

.980000 

.539511 

. 1223B7E+07 

.980000 

3 . 60579 

710459. 

.980000 

. 560402 

. 129311E+Q7 

.980000 

3.74542 

698347 . 

.980000 

.582103 

. 136104E+07 

.980000 

3 . 89045 

686362 

.980000 

. 604644 

142713E+07 

.980000 

4.04111 

674540 . 

.980000 

.628058 

. 149084E+07 

.980000 

4.19759 

662918. 

.980000 

. 652379 

. 1S5160E+D7 

.980000 

4.36014 

651S31 

.980000 

. 677641 

. 160B8BE+07 

.980000 

4 . 52898 

640408. 

.980000 

.703882 

. 166214E+07 

.980000 

4.70436 

629579 . 

.980000 

.731139 

. 171085E+07 

.980000 

4 . 88653 

619069 

.980000 

. 7S9451 

. 175451E+07 

.980000 

S . 07575 

608902. 

.980000 

. 788860 

. 179265E+07 

. 9BOOOO 

5 . 27230 

599099 . 

.980000 

.819408 

1B24B6E+07 

.980000 

5 . 47646 

589680 . 

.980000 

.851138 

. 1B507SE+07 

.980000 

5 . 6BBS3 

S80661 . 

.980000 

.884097 

. 1B7002E+07 

.980000 

5.90881 

572058 

.980000 

.918333 

. 188241E*D7 

.980000 

6.13762 

563882 . 

.980000 

. 953894 

. 1BB774E+07 

.980000 

6.37529 

556146 

.980000 

.990832 

. 188592E407 

.980000 

6.62217 

548860. 

.980000 

1 . 02920 

. 187692E+07 

.980000 

6.87860 

542032. 

.980000 

1.06905 

. 186080E+07 

.980000 

7.14497 

535671 . 

.980000 

1 . 11045 

. 183772E+07 

.980000 

7.42164 

529779 . 

.980000 

1 . 15345 

. 180791E+07 

.980000 

7.70904 

524366 . 

.980000 

1  . 19812 

. 177165E+07 

.980000 

B. 00756 

519434. 

.980000 

1.24451 

.  172934E+07 

.980000 

8.31764 

514987. 

.980000 

1  29271 

. 16B142E+07 

.980000 

8 . 63973 

511027. 

.980000 

1  34277 

. 162840E+07 

.980000 

8 . 97429 

507559 

.980000 

1  39476 

1S7QB3E*07 

.980000 

9.32181 

504583 

.980000 

1 . 44877 

. 1S0931E+07 

.980000 

9 . 68279 

502102. 

.980000 

1  50487 

.  144446E+07 

.980000 

10.0577 

500117. 

.980000 

1.56315 

. 137693E+07 

.980000 

10.4472 

49B631 . 

.980000 

1 . 62368 

. 130737E+07 

.980000 

10.8518 

497644 . 

9BOOOO 

1  68655 

. 123642E+07 

.980000 

11.2720 

497158. 

.980000 

1.75186 

.  U6472E+07 

.980000 

11.7085 

497174. 

.980000 

1.81970 

. 109286E+07 

.980000 

12.1619 

497695. 

.980000 

1 .89017 

. 1D2142E+07 

.980000 

12.6328 

498720 . 

.980000 

1 . 96336 

950925 . 

.980000 

13.1220 

500252. 

.980000 

2 . 03939 

881865. 

.980000 

13.6301 

502295 . 

.980000 

2.11836 

866973 . 

.980000 

14.1579 

S04848 . 

.980000 

2.20039 

858143. 

.980000 

14.7062 

507912. 

.980000 

2.28560 

848746. 

.980000 

15.2757 

511491 . 

.980000 

2.37411 

838833 . 

.980000 

15.8672 

515586 . 

.980000 

2.46604 

828454. 

.980000 

16.4816 

520200 

.980000 

2 . 56153 

817639. 

.980000 

17.1199 

525333 . 

.980000 

2.66073 

806500. 

.980000 

17.7828 

530989 . 

.980000 

2 . 76376 

795035. 

.980000 

18.4714 

537170 . 

.980000 

2.87078 

783312. 

.980000 

19 . 1867 

543874 . 

.980000 

2  98195 

771385. 

.980000 

19.9297 

551105 . 

.980000 

3.09742 

7S9305 . 

.980000 

20.7014 

558863 

.980000 

3  21736 

747125. 

.980000 

21.5031 

567146. 

.980000 

3.34195 

734891 . 

.980000 

22.3357 

575956 . 

.980000 

Laabda 

JUnga-Bandvidth 
Product  (tt-Hs) 

Paak 

Laabda 

Rang*- Bandwidth 
Product  (ft-Hs) 

Paak 

23.2007 

585293 . 

.980000 

155.060 

. 147214E+07 

.980000 

24.0991 

59S155. 

.980000 

161.065 

. 149S96E+07 

.980000 

25 . 0323 

605537. 

.980000 

167.302 

. 152016E+07 

.980000 

26.0016 

616441 . 

.980000 

173.780 

1S447SE+07 

.980000 

27.0085 

627857 . 

.980000 

180.510 

. 1S6974E*07 

.980000 

2B . 0544 

639780 . 

.980000 

187.500 

. 1S9514E*07 

.980000 

29.1407 

652212 . 

.980000 

194.760 

. 162094E+07 

.980000 

30.2692 

665135 

.980000 

202.302 

. 164717E*07 

.980000 

31.4413 

678555 

.980000 

210.136 

. 167382E+07 

.980000 

32.6S88 

692445 . 

.980000 

218  273 

. 1700B9E+07 

.980000 

33  9235 

706804. 

.980000 

226  725 

. 172841E+07 

.980000 

35.2371 

721629. 

.980000 

235.505 

. 175638E+07 

980000 

36.6016 

736895 . 

.980000 

244.625 

. 178479E+07 

.980000 

38.0190 

752S82 . 

.980000 

254.097 

. 181366E+07 

.980000 

39.4912 

768686 . 

.980000 

263.937 

. 1B4300E+07 

.980000 

41  0205 

785199. 

.980000 

274.158 

.  1B72B2E+07 

.980000 

42.6089 

802083. 

.980000 

284.774 

. 190312E+07 

.980000 

44  2589 

819329 

.980000 

295.802 

. 193391E+07 

.980000 

45.9728 

836941 . 

.980000 

307.256 

. 196519E*07 

.980000 

47.7S30 

854869 . 

.980000 

319.154 

. 199699E*07 

.980000 

49.6022 

873111 . 

.980000 

331.513 

. 202929E*07 

.980000 

51.5229 

891650 . 

.980000 

344.350 

. 206212E+07 

.980000 

53.5181 

910458. 

.980000 

357 . 685 

. 209548E*07 

.980000 

55.5905 

929514. 

.980000 

371.535 

.212939E+07 

.980000 

57.7431 

948860 . 

.980000 

385 . 923 

. 216383E+07 

.980000 

59.9791 

968364 . 

.980000 

400.867 

. 219B84E+07 

.980000 

62.3018 

988119 

.980000 

416.390 

. 223441E+07 

.980000 

64  7144 

. 100803E+07 

.980000 

432.514 

. 227056E+07 

.980000 

67  2203 

. 102820E+07 

.980000 

449.262 

. 230730E+07 

.980000 

69  8233 

104849E+07 

.980000 

466.659 

.  234462E+0? 

.980000 

72.S271 

106796E*07 

.980000 

4B4.731 

. 238255E+07 

.980000 

75  3357 

106524E+07 

.980000 

503.501 

.  2421 09E+07 

.980000 

78  2529 

. 110279E*07 

.980000 

522 . 998 

. 246026E+07 

.980000 

81  2831 

.  112064E+07 

.980000 

543.250 

.  250007E+07 

.980000 

84.4307 

. 113877E+07 

.980000 

564  287 

. 254051E+07 

.980000 

87.7002 

.  115719E+07 

.980000 

586 . 138 

. 258161E+07 

.980000 

91 . 0963 

.  117591E+07 

.980000 

608 . 835 

. 262337E+07 

.980000 

94  6238 

.  119493E*07 

.980000 

632.412 

. 2665B2E+07 

.980000 

98.2880 

. 121426E*07 

.980000 

656 . 902 

. 270895E+07 

.980000 

102.094 

. 123391E+07 

.980000 

682 . 339 

.27S277E+07 

.980000 

106.048 

. 1253B7E+07 

.980000 

708.761 

. 279730E+07 

.980000 

110. 1S4 

. 127416E+07 

.980000 

736.208 

.284256E+07 

.980000 

114  420 

. 129477E+07 

.980000 

764.717 

. 288855E+07 

.980000 

118.850 

.  131S71E+07 

.980000 

794.330 

. 29352BE+07 

.980000 

123  453 

. 133700E+07 

.980000 

82S.089 

298277E*07 

.980000 

128.233 

. 135B63E+07 

.980000 

857 . 039 

.  303102E+07 

.980000 

133.199 

. 13B061E*07 

.980000 

890 . 227 

. 3080065*07 

.980000 

138  357 

.  140295E+07 

.980000 

924.699 

.312989E+07 

.980000 

143.714 

.  142564E*07 

.980000 

960.508 

. 318D52E+07 

.980000 

149.280 

144871E*07 

.980000 

997.702 

323197E+07 

.980000 

Lambda 

Ranga-Bandvidth 
Product  (ft-H*) 

Paak 

'-ambda 

Ranga-Bandvidth 
Product  (ft-Hz) 

Paak 

.519397 

. 265469E+07 

.990000 

3 . 47136 

.  166266E+07 

.990000 

.539511 

. 2815B3E+07 

.990000 

3 . 60S79 

. 163460E+07 

.990000 

.560402 

. 297515E+07 

.990000 

3.74542 

. 160673E+07 

.990000 

.582103 

313144E+07 

.990000 

3.89045 

. 15791SE+07 

.990000 

.604644 

.328349E+07 

.990000 

4.04111 

. 1S519SE+07 

.990000 

. 628058 

. 343006E+07 

.990000 

4.19759 

. 152522E+07 

.990000 

. 652379 

. 356986E+07 

.990000 

4.36014 

. 149902E+07 

.990000 

.677641 

. 370166E+07 

.990000 

4 . 52898 

. 147343E+07 

.990000 

.703882 

. 3B2419E+07 

.990000 

4.70436 

.144BS1E+07 

.990000 

.731139 

. 393625E+07 

.990000 

4 . 88653 

. 142433E+07 

.990000 

.759451 

. 403671E+07 

.990000 

5 . 07575 

. 140094E+07 

.990000 

.788860 

. 412447E+07 

.990000 

5.27230 

. 137838E+07 

.990000 

.819408 

. 4198S7E+07 

.990000 

5 . 47646 

.  135671E+07 

.990000 

.851138 

425814E+07 

.990000 

5 . 68853 

. 133596E+07 

.990000 

.884097 

. 430246E+07 

.990000 

5.90881 

. 131617E+0? 

.990000 

. 918333 

.433097E+07 

.990000 

6.13762 

. 129736E+07 

.990000 

.953894 

. 434323E+07 

.990000 

6 . 37529 

. 1279S6E+07 

.990000 

.990832 

. 433904E+07 

.990000 

6  62217 

. 126280E+07 

.990000 

1.02920 

. 431B34E+07 

.990000 

6.87860 

. 124709E+07 

.990000 

1 . 06905 

. 428126E+07 

.990000 

7.14497 

. 12324SE+07 

.990000 

1.11045 

. 422B17E+07 

.990000 

7.42164 

. 121B90E+07 

.990000 

1 . 15345 

. 415956E+07 

.990000 

7.70904 

. 120644E+07 

.990000 

1 . 19812 

. 407613E+07 

.990000 

8.00756 

.  119509E+07 

.990000 

1 .24451 

. 397B79E+07 

.990000 

8.31764 

. 1 18486E+07 

.990000 

1.29271 

. 386854E+07 

.990000 

8 . 63973 

.  117575E+07 

.990000 

1 . 34277 

. 374655E+07 

.990000 

8.97429 

. 116777E+07 

.990000 

1 . 39476 

. 361409E+07 

.990000 

9  32181 

116092Ef07 

.990000 

1  .44877 

347256E+07 

.990000 

9 . 68279 

.  115522E+07 

.990000 

1 .50487 

. 332335E+07 

.990000 

1 0 . 0577 

.11S06SE+07 

.990000 

1.56315 

316799E+07 

.990000 

10.4472 

. 114723E+07 

.990000 

1 . 62368 

. 300795E+07 

.990000 

10.B518 

. 114496E+07 

.990000 

1 . 68655 

284472E+07 

.990000 

11  2720 

. 114384E+07 

.990000 

1.75186 

,267974E*07 

.990000 

11  7085 

.  1143B8E+07 

.990000 

1 .81970 

.251440E+07 

.990000 

12.1619 

.  114508E+07 

.990000 

1.89017 

. 235003E+07 

.990001 

12.6328 

. 114744E+07 

.990000 

1 . 96336 

21878SE+07 

.990000 

13.1220 

. 115096E+07 

.990001 

2.03939 

. 202B96E+07 

.990000 

13.6301 

.  115566E+07 

.990000 

2.11836 

. 199470E+07 

.990000 

14.1579 

. U6153E+07 

.990000 

2.20039 

. 197438E^B7 

.990000 

14.7062 

.  116858E  +  07 

.990000 

2.28560 

. 19S276E+07 

.990000 

15 . 2757 

.  1176B2E+07 

.990000 

2  37411 

. 192995E+07 

.990000 

15.8672 

. 118624E+07 

.990000 

2 . 46604 

. 190607E+07 

.990000 

16.4816 

. 1 19686E+07 

.990000 

2  561 S3 

. 188124E+07 

.990000 

17.1199 

.  120B67E+07 

.990000 

2  66073 

. 1B5556E+07 

.990000 

17.782C 

. 12216BE+07 

.990001 

2 . 76376 

. 182919E+07 

.990000 

18.4714 

.  123590E+07 

.990001 

2  87078 

. 180221E+07 

.990000 

19.1867 

. 125132E+07 

.990000 

2.98195 

. 177477E+07 

.990000 

19.9297 

. 126796E+07 

.990000 

3 . 09742 

. 174698E+ 07 

.990000 

20.7014 

. 12B581E+07 

.990000 

3.21736 

. 171896E*07 

.990000 

21.5031 

■130487E+07 

.990000 

3  34195 

169081E+07 

.990000 

22.3357 

. 132514E+07 

.990000 

Laabda 

Range-Bandwidth 
Product  (ft-Hz) 

Peak 

bttbda 

Range- Bandwidth 
Product  (ft-Hz) 

Peak 

23.2007 

. 134662E+07 

.990000 

155.060 

. 338705E+07 

.990000 

24.  om 

. 136931E+07 

.990000 

161.065 

. 344184E+07 

.990000 

25.0323 

139320E+07 

.990000 

167.302 

. 349753E+07 

.990000 

26.0016 

. 14182BE+07 

.990000 

173.780 

3SS411E+07 

.990000 

27.0085 

. 144455E+07 

.990000 

180.510 

. 361161E+07 

.990000 

28 . 0544 

.  147198E+07 

.990000 

187.500 

. 367004E+07 

.990000 

29.1407 

. 150058E+07 

.990000 

194.760 

. 372940E+07 

.990000 

30.2692 

.1S3032E+07 

.990000 

202.302 

. 378974E+07 

.990000 

31.4413 

. 156119E+07 

.990000 

210.136 

. 38510SE+07 

.990000 

32 . 6588 

•15931SE+07 

.990000 

218.273 

. 391335E+07 

.990000 

33  9235 

. 162619E+07 

.990001 

226 . 725 

397666E+07 

.990000 

35.2371 

.  166030E+07 

.990000 

235.505 

404101E+07 

.990000 

36.6016 

. 169542E+07 

.990000 

244  625 

410637E+07 

.990000 

38.0190 

. 173151E+07 

.990000 

254.097 

. 417280E+07 

.990000 

39.4912 

. 176B56E+07 

.990000 

263 . 937 

. 424031E+07 

.990001 

41.0205 

. 1B0656E+07 

.990000 

274.158 

. 430B92E+07 

.990000 

42 . 6089 

. 184540E+07 

.990000 

284.774 

. 437863E+07 

.990000 

44.2589 

.  188508E+07 

.990000 

295.802 

444947E+07 

.990000 

45 . 9728 

. 192S60E+07 

.990000 

307.256 

. 4S2144E+07 

.990000 

47.7530 

. 196685E+07 

.990000 

319.154 

. 459459E+07 

.990000 

49.6022 

200BB2E+07 

.990000 

331  513 

. 466892E+07 

.990000 

SI .5229 

20S147E+07 

.990000 

344.350 

. 474446E+07 

.990000 

53.5181 

. 209475E+07 

.990000 

357 . 685 

. 482121E*07 

.990000 

55.5905 

.213B59E+07 

.990000 

371. S3S 

. 4B9921E+07 

.990000 

57  7431 

. 218310E+07 

.990000 

385  923 

497846E+07 

.990000 

59  9791 

. 22279BE+07 

.990000 

400.867 

.505900E+07 

.990000 

62.3018 

. 227343E+07 

.990000 

416.390 

. 514085E+07 

.990000 

64  7144 

.231924E+07 

.990000 

432.514 

S22402E+07 

.990000 

67.2203 

. 236565E+07 

.990000 

449 . 262 

. 530B54E+07 

.990000 

69.8233 

.241232E+07 

.990001 

466.659 

539441E+07 

.990000 

72.5271 

. H4S712E+07 

.990000 

484.731 

. 54B168E*07 

.990000 

75.33S7 

.249688E+07 

.990000 

503.501 

557035E+07 

.990000 

78 . 2529 

. 253727E*07 

.990000 

522  998 

. 56604BE+07 

.990000 

81.2831 

.257832E+07 

.990000 

543.250 

. 575205E+07 

.990000 

84.4307 

. 262003E+07 

.990000 

564 . 287 

584510E+07 

.990000 

87.7002 

. 266242E+D7 

.990000 

SB6.138 

593967E+07 

.990000 

91.0963 

. 270549E*07 

.990000 

608.835 

603575E+07 

.990000 

94  6238 

.274925E+07 

.990000 

632  412 

. 613341E+07 

.990000 

98  2880 

. 279373E+07 

.990000 

656.902 

. 623264E+07 

.990000 

102.094 

. 283893E+07 

.990000 

682.339 

.633347E+07 

.990000 

106.048 

. 2B84B6E+07 

.990000 

708.761 

. 643593E+07 

.990000 

110.154 

. 293153E+07 

.990000 

736.208 

654005E+07 

.990000 

114.420 

. 297896E+07 

.990000 

764.717 

. 6645B7E+07 

.990000 

118  850 

.  302714E+07 

.990000 

794.330 

675338E+07 

.990000 

123  453 

. 307612E+07 

.990000 

825.089 

. 686264E+07 

990000 

128  £33 

.  312589E+07 

.990000 

857.039 

. 697365E*07 

.990000 

133  199 

317645E407 

.990000 

890 . 227 

. 70864BE+07 

.990000 

138.357 

322784E+07 

.990000 

924  699 

. 720113E+07 

.990000 

143.714 

. 328006E+07 

.990000 

960.508 

731761E407 

.990000 

149.280 

. 333313E+07 

.990000 

997.702 

. 743600E+07 

.990000 
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